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ABSTRACT 


In this dissertation an automatic economic design 
of linear discrete structures is presented. Two types of such 
structures are considered: (a) pin- jointed plane and space 
trusses, and (h) rigid jointed plane frames which include 
continuous heams, gable bents, multi-storeyed buildings and 
arches. For both types of structures, the analysis is done 
by point iteration of the generalised slope deflection 
equations. 

For both cases feasibility of stmctural design 
is examined. Pin-jointed trusses are designed for single as 
well as multi-load conditions by an iterative procedure such 
that the volxme of the structural members is a minimum. The 
resulting structure need not be a fully stressed one. 

Rigid jointed structural members are designed 
to achieve the minimum sectional area by iterating on the 
depth of the section. The sections used are prismatic 
I-section members. 

The design procedure is completely automated 
once the geometiy and design loads are specified. 



The structure is first analysed by approximate 


method and approximate section properties are arrived at. 

The properties are, in turn, utilized in the more exact 
analysis and the design cycle is repeated until an economic 
depth and corresponding section property for eveiy monber 
is obtained. 

Several problems are illustrated and discussions 
of the results are presented. 



NOMENCLATTJEE* 
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BS 

CM 

CM 
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ij 
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I_ 
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Area of cross-section 
Area of flange 

Spread of flange 
Interaction value 
Beam bending moment 
Beam shear force 
Column moment 

Colximn bending moment for rth storey 

Depth of section 
Young’s modulus 

Allowable stress in axial and bending 

Actual stress due to axial force and bending 
moment 

Force in I, J member 
Depth of a section 

Moment of inertia and used as a subscript 
Moment of inertia of the columns of rth storey 

Column moment of inertia 

Cirder moment of inertia 

Moment of inertia about k axis 


J Polar moment of inertia 

*Sso see nomenclature for computer programmes 
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K 




L 



Stiffiiess matrix 
Element stiffness matrices 

Span of member 

Coliamn and girder lengths 

Bending moment 

Moment of resistance of 42^'th storey 


N 

P 



Column noimal force for rth storey 

Actual force 
Load at a joint 

Distance from C-G. of a column of ith storey 
Reaction of ij joint 
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or 



al 


Critical stress to prevent failure due to 
buckling 

Critical bending stress to prevent failure 
due to bending buckling 
Transformation matrix 

Thickness to prevent bucklixig due to shear 

axial and bending moments 

Jont displacement 

Shear force 

Shear at rth storey 

Critical value of shear stress to prevent 
buckling 

Allowable shear stress 
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max 

Design shear force 

V 

Inspan load 

Y. 
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Distance of a section fron neutral axis 

Z 

Section modulus 

09 0 • • 9 0 

’ r 0 

Rotation of a joint 


Displacement of a joint 


Poissions Ratio 


Storey sway angle 

cr 

Stress 



CHAPTER I 


IMTRODUCTION 

The structural design may he defined as the process 
of determining the most suitable geometry of a struottire to 
satisfiy a set of purposes and to dimension the structural ele- 
ments of which it is composed of. A most suitable geomet3:y 
has to satisfy the functional, economic, social, aesthetic 
and many more requirements. Once the configuration is obtained, 
the forces or loads that the structure has to resist is to be 
determined. These loads are tabulated in the Codes of Speci- 
fications and does not require much attention for the conven- 
tional type of structures, but may really be a difficult 
problem for special and rare type of structures like space 
vehicles, supersonic jet planes etc. It is beyond the scope 
of this thesis either to design the geometry or to determine 
the design forces, instead both configuration and design load 
systems are assumed. So the problem dealt here is to propo3^“ 
tion the elements of the structure subjected to a given set 
of loads. The structures considered are the rigid jointed 
plane frames, like multi-storeyed building, gable frames and 
beams? and pin jointed plane and space trusses. Plane arches 
may approximately be solved if the curved geometry of the arch 
is replaced by an equivalent linear gecmetry. The length of 
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the arch along the centre line may he divided into a set of 
linear elements for this purpose. 

1.1 Design of Elements 

There are two distinct phases in the design of 

I 

the meanhers of a structure: 

Phase 1; 

The stress and deformation analysis of the entire 
structure has to he completed and design forces for elements 
obtained. The stress analysis depends upon a set of equations. 
Henceforth these equations will he called as the equations of 
analysis. These are the equilibrium equations which state that 
the internal forces in the structure and the external loads 
must he in equilibrium and the equations of oompat ability which 
states that the deformed msnhers of the structure must continue 
to fit together, so that the deformations are geometrically 
compatible. There is another set of equations called the 
stress strain laws. Since the structures dealt here are linear 
structures these equations are automatically satisfied and do 
not appear in the analysis. 

Phase 2: 

The members are to be designed for ihe above forces 
of the analysis such that the stress at the hi^ly stressed 
sections should be within the allowable limits. The deforma- 
tions at any point along the span of the member should also 
be within the permissible limits. This is the elastic design 
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of structures. For future references these ■will he known as 
the constraints of •the design. Several more constraints on 
depth, width and thickness of the section may he laid do'wn. 

In short any valid design of elsaents must 
satisfy the equations of analysis and the constraints of 
design. 

The above design does not ensure that the volume 
of the structure is a minimum. If in addition to the above 
equations and constraints one more constraint of the volume 
i.e. the total volume of the structure he a minimum is laid 
down the problem becomes a minimum weight design problem. 
Mathematical programming has to be resorted to for the solu- 
tion of thisproblem. If all the equations of analysis and 
constraints of design are linear the» it is a linear pro- 
gramming problem otherwise it is a non-linear one if the 
equations or constraints are not linear. Elastic design is 
a non-linear programming problem because the equations of 
equilibrium are non-linear. These equations for trusses has 
the form 

f(€r*. A.) » P. .. 1*1 

JU JL X 

whe3?e both the member stress area of the member 

'Aj^' are unknowns and they occur in pairs making it a non- 
linear set of equations. is the external load vector. 

This lype of design where member parameters and manber 
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stresses are detemined to minimize the volume (also called 
the objective f\aiction in ihe literature) is the direct design 
of structures. 

The above design is based on the stress criterion. 
There is another type of design in which the basis is the ulti- 
mate strength of the structure. This approach goes under the 
name limit design or -tlastio design of stiructures. Like a 
direct elastic design there is also a direct plastic design of 
structures. The plastic design of structures comes under 
linear programming. 

Althou^ the direct design is an excellent method 
the labour involved is extremely high and with the present com- 
putation facility only small structures like a two storey portal 
frame, gables, small trusses with only a few panels are attempted. 

Since contemporary structures are highly compli- 
cated and are unusually large like a hundred storey building or 
a long truss bridge, the above direct desior with the present 
knowledge and computation facility becomes highly involved as 
the number of equations of analysis and constraints of design 
are unusually large. Much work is being currently carried cut in 
direct design of structures. F. Moses (l)* has transformed the 
non-linear elastic design of a truss to a step by step linear 


*Figui*e in brackets refe3?s to reference number in 


the 


Bibliography. 
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design problem by Taylor se3rf.es expansion. L.C. Schmidt (2) 
has developed methods for the direot elastic design changing 
even the geometsy of the st 33 ucture. Bigelow (j) has woiiced on 
the linear p3rog3?ainming pixjblem applied to plastic design. 

Other pioneers on this field are Prager (4)» Brown (5)» Save (6) 
and Koopman (7). 

In contrast to the direct design there is also the 
indirect design also called the interative design of st3nictures. 
In this method first the analysis is completed and the elements 
are propoortioned acoo33ding to the design constraints. The cycle 
of analysis and design are 3:epeated until the sections p3:opea?- 
ties remain constant i.e. the iteration on design conveiiges. 

Malcolm (s) puts the above iteirative design in a 
sequential fozm ast 

(1) Estimation (i.e. guessing) of the member sizes. 

( 2 ) Analysis for stress and defomation for the 
above proposrbions. 

( 5 ) Checking of the sections for strength stability 
and seorviceability for the forces of the above 
analysis. 

( 4 ) Successive modifications of the 03rf.ginal guess 

for member sizes until convergence in the member 
dimension is achieved. 

For a clear mderstanding of the porilnciples undej?- 
lining the four steps of an iterative design as suggested by 
Malcolm the second step is treated first. 



1.2 Stress and Deformation Jinalysis of a Structure 

This is the most important and involved procedure 
in the iterative design of structures. This is concerned with 
the evaluation of forces and deformation at every point of the 
stiruoture. 

There ane two fundamental methods of analysis, the 
method in which the laws of equilihrium are applied first leads 
to the equations of geometric compatibility and is known as the 
Compatibililgr or the Force or the Flexibility method and the 
method in which the conditions of compatibilities are used first 
gives rise to the equations of joint equilibrium and is called the 
Equilibrium or the Displacements or the Stiffness method. 

Both the methods finally lead to a set of simulta- 
neous equations which may be the compatibility or the Equilibrium 
equations. 

Although any of the methods can be used for the 
pnalysis of any structur, it is obvious that for a given struc- 
ture the method which is less involved should be chosen. 

1 . 5 Force Method 

For hand calculation ELexibilily method is less 
involved and is more widely used for structures with low degree 
of statical indeterminacy, as the set of equations to be solved 
are less in number. The persons associated with this method are 
Argyris (9) and P.B» Morioe ( 10 ) althovigh the method was first 
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formulated in 1888. The method of consistent deformations 
and all Energy approaches (11) of analysing structures come * 
under this method. Basically this method is the seme as the 
method of Mohr (12) who formulated the basic equations and 
called them as the simultaneous condition equations of an in- 
determinate structure which egresses the geometric coherence 
of the structure (i.e. the geometric compatibility of the 
structure) . 

Up to the present time Force method was popular 
but in the machine computation this method is not generally 
used, althou^ for trusses the method yield very few equations, 
because the process of generating these equations in a computer 
involve two major difficulties: 

(1 ) In the Ibrce method to generate the equations a 
stable released structure is necessary which is to be obtained 
out of the real structure by removing all the redundant s. For 
a oompler structure the determination of the degress of redun- 
dancies and a choice of a released structure which may save a 
lot of computation time is extremely difficult as this is a 
structure dependent phenomenon and varies from structmre to 
structure. 

( 2 ) The compatibility equations are the equations of 
deformations. Even when all the member forces are known to 
calculate the defoliations and hence to generate each ccmpati- 
bility equation of the Ibrce method all the members of the 
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structurehr.ve to be scanned. For hand calculations many short- 
cuts may be resorted to, to simplify the above generation of 
the final equations depending upon the structure. Since a 
computer programme is generally written for all types of struc- 
tures such labour and computation saving is not possible 
which make the computation tedious and time consuming. Th is 
difficulty for future reference will be known as the struc- 
ture wide effect of Ebrce method. 

In short, Force method is generally not recommended 
because of the difficulties involved for the generation of the 
equations, although the number of these equations for low in- 
determinate structures are very less. For example, for a tiruss 
with two redundants and 50 joints such equations are only two 
whereas for displacement method the number may be 100 or 150 
according to the dimensions (plane or space) of the truss. 

1.4 Displacement Method 

The Displacement method of structural analysis 
is basically due to Maney (l 5) who published the slope deflec- 
tion equations in 1915* The number of equations grow rapidly 
as the structures gets complicated, for example a four storey 
and four bay plane building frame considered to be a very 
fflnall structure, gives rise to 52 simultaneous equations consi- 
dering only two displacements (joint rotation and joint trans- 
lations) for each of its I6 joints. It is extremely difficult 
to solve the 52 equations by hand. Hence before the intro- 
duction of calculating machines various relaxation and 



iterative procedures are developed to solve the structures by - 
hand. Few of the notable methods for analysing structures by 
hand are given below. 

1 . Moment Distribution 

This is a relaxation method and is due to Hardy 
Gross who published this in 1950* this method the build- 
ing frame is held against sway at each floor level and joint 
equilibriums are satisfied by distributing the mbalanced 
moments among all the members meeting at the joint according 
to their stiffnesses. Storeys are then allowed to sway. and 
member end momentsdetermined for each sway. These moments are 
balanced and storey shear equations wiitten for each storey 
which gives the sway. By superposing the solution of unsway 
and sway cases the final moments are calculated. This method 
is unsuitable for multi-storeyed frames as the number of moment 
distributions increases according to the niamber of storeys but 
are powerful for non-translatory structures. 

2. Kanis Method ( 14 ) 

6. Kani developed an iteration method to solve 
the slope-deflection equations. In this method both non-sway 
equations and storey shear equations are balanced together. 

This is an excellent contribution but limited to building 
frames with constant storey height for a storey. 

3 . Takabeya Method ( 15 ) 

Also known as the Japanese method iterates 



the slope deflection equations at both far end and near end 
simultaneously for a few cycles and then distributes the un- 
balanced moment at a particular joint according to the prin- 
ciple of moment distribution. This seoas to be a powerful 
method as Takabeya has solved a two hundred storey and twenty- 
seven bay building frame by hand within seventy-eight hours. 

The ccmplexity of modem structures and the use 
of digital computer has discarded the above methods. The 
equations of the Displacement method are written in the matrix 
form and is the most popular method in machine computation. 

It does not have any of the difficulties encountered in the 
Force method. Each equation of this method deal only with 
the number of members meeting at a joint and ere are comple- 
tely localized to the joint and the for end joints of the 
member meeting at the joint. This gives rise to a regular 
band for ''most of the symmetric and regular structures. Since 
redundants of the structures are not needed at any phase of 
the analysis the difficulty of detennining a basic structure 
necessary in the Force method is avoided. The popularity 
of this method is due to the above two advantages, althou^ • 
the final equations of ULsplacement method for trusses are 
sometimes unusually large than the equations of the porce 
method. 

There are two ways of ainriving at the final 
equations of the Displacement method. The first method is 
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based on the principle of virtual woik and involves the 
multiplication cf a number of matrices the dimensions of which 
are much larger (7 to 8 times larger) than the dimension of 
the final equations and is seldom used in the computers. 

Semih S. Tezoan (l6) and Pestal (17) use this method to arrive 
at the final equations. Semih S. Tezcan (l6) assumes the sway 
of a storey constant i. e. the axial deformations are neglected 
where as Pestel (l7) includes axial deformations of beams and 
coltimns. In the other method joint equations are directly 
written and is known as the method of merging. 

In 1962 Eiseman and Namyet (18) published the 
above approach for space frames. They have generalized the 
slope-deflection equations for discrete space structures. 

They have developed a computer programme called ¥RM 
((^Rone ANalysis) in which the total structure is divided into 
several units and eliminating some of the joint displacements 
while proceeding from one unit to other they are able to solve 
huge structures in computers using external storage. 

Semih S. Tezoan (20) has also used the generalised 
slope deflection equations to arrive at the final equations. in 
matrix form. Then breaking up the structure stiffhess matrix* 
into two trangular matrices and using Cholesky's square root 
scheme, has shown that without the help of external storage, 
1000 equations with a band width of 27 can be solved by a 
IBM 7090 computer having 52000 word capacity which othersise 
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can invert only a I50 x I50 matrix. With the help of external 
storage he has handled as hi^ as 3500 equations. An another 
paper Tezcan (21 ) has fomulated the joint equilibrium equa- 
tions using slope objlection equations andneglecting the axial 
deformations. Here he has formulated the joint and storey 
equations separately. 

E.W. Clough and Ian. P. King (22) have developed 
a programme in which the floor displacements are grouped to- 
gether, thereby arriving at a three band stiffness matrix for 
each frame of the building. The effect of shear walls in 
between two girders are taken in formulating the slope deflec- 
tion equations and they have solved the final matrix by parti- 
tioning technique. 

J.E. Goldberg (23) has contributed a method to 
calculate stresses due to lateral loads in multi-storeyed 
building frames. He has modified the slope deflection equa- 
tions to take into account shear stiffners in the form of walls 
or sha.ft8« The equations are generated from top floor of the 
building and while proceeding to the bottom most floor, the 
displacements of the top floors are eliminated and where the 
bottom most floor is reached the equations contain only the 
displacements of this floor. This displacements are deter- 
mined and then the building is scanned from bottom to top 
to yield the total displacement matrix. 


Other people associated with this method are 
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Archer (24) » Benscotter (25), Behke (26), Klein (27)5 Martin 
(28), Turner (29) ad et. 

'basically most of the works deal with two 
fundamental procedures: 

1. The Model: 

Here the structure is idealised namely - shear 
wall are introduced, building is visualized as a space, tier 
or a plane frame. Trusses are introduced to repi?esent the 
slabs and walls and all types of improvement is thou^out of 
to represent the real structure by a model which may res^ble 
the r.ctufil stiuoturo closely. 

2. The Computation: 

The above model leads to a set of equations and 
the researcher developed means of accommodating and solving 
these equations with the present digital computers having a 
limited memory. 

Most of the modem structures are r&ry complex 
in nature and which need not give a banded set of equations 
like rings and space trusses. To solve these structures the, 
above procedure is sometime quite involved with the present 
computation facility. For such structure an iterative tech- 
nique is developed. In this procedure joint stiffness matrix 
is generated and sisplacements of the joint is calciilated. 

The procedure is repeated for all joints several times until 



the displacements have a stationary value. This type of 
solving the structures saves the memory of a computer as 
the only storage necessary is the joint stiffness matrix. 

To accelerate the ccnvergancy initial values of joint dis- 
placements are calculated hy approximate methods. This is 
similar to Gauss 5e idel method in -which only the required 
number of equations necessary are generated and stored. This 
method is ajplicahle -to any type of ir3?egular structure and 
there is no need for a banded stiffness matrix. 

Once the analysis is complete next step is the 
design of a struc-ture. As already stated the most important 
factor of a design is the economy. Since the structures dealt 
are rigid jointed and pin- join ted structures and eanh need 
different treatment they are taken up separately. 

1.5 Pin- jointed Structures 

In 1900 Gilley (50) pointed out that for an in- 
determinate truss under one load condition, the minimum -weight 
structure is a statically determinate structure in which the 
cross-section dimension of each member is selected to achieve 
the allowable stress for that member. For multiple load con- 
dition he used the concept of a fully-stressed design. A 
fully-stressed design is defined as one in which each msaber 
reaches its allowable stress under at least one of the load 
system. In 1954 Sved (31) also proved that under one load 
condition the minimum weight truss is a determinate one. 



15 


In 1958 Schmidt (52) woriced on the indeter- 
minate truss and noted that there are several fully-stressed 
designs for a truss under the load conditions and the least 
weight fully-stressed structure is the minimum weight struc- 
ture. But in 1960 L.A. Schmit ( 35 ) solved the above problem 
using mathematical programming and came to the conclusion that 
minimum weight structures are not necessarily fully-stressed 
designs. 

In 1964 Eazani (34) examined the relationship 
of the minimm weight design and the fully-stressed design 
and concluded that fully-stressed design need not converge to 
a minimum weight design. In 1986 Filliam Weaver ( 35 ) developed 
an automatic design procedure for space trusses. This is a 
fully-stressed design and also includes dynamic effects like 
the resonant frequency requirements. 

Most of the above structures have a fixed geometry. 
Schmit extended the concept of overall geometric configuration 
as a design parameter. A direct design, to evalxxate the geo- 
metric configuration to carry the load was taken up by Michell 
( 36 ) in 1904. His structures are known as Michell structures, 
the geometiy of these structures are such that th^ are not 
practically usable. But the theorems and procedure of Michell 
is extensively used in the aircraft industry. 

In this thesis minimum wei^t structures are 
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designed by an iterative procedure. (The concept of a fully- 
stressed design for one load condition is explained in the 
Chapter 'Feasibility of a Structural Design'. It is proved 
that for an indeterminate structure with one load condition 
a full-stress design is not at all possible as it violates 
the fundamental compatibility equations of a pin-jointed 
structure. It is also shown that even for multiple load con- 
dition, a fully stressed design may not be possible. The 
design procedure for these structures are dealt- in this Chapter 
and in the section 4*15 of Chapter I?. 

1.6 Rigid- jointed Structures 

If the research on the analysis of such struc- 
tures are compared with the research on design it can be 
clearly seen that analysis rather than design has consumed 
far greater share of attention althou^ analysis is only a 
part of the design procedure. Sophisticated computer pro- 
grammes are developed for the analysis of discrete and conti- 
nuous structures but for design one depends upon the age old 
methods of cut and try procedure. 

In 1959 Grinter (5?) pointed out this fact and 
introduced a concept of automatio desiga. of beams and frames. 

In this procedure he tried to distribute the fixed end modulus 
instead of the fixed end moments in the moment distribution 
procedure. At intermediate steps he charged the distribution 
factors so that at the eiad of the distribution the balanced 
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modulus and the actual modulus agreedwithin a short margin. 

He has also shown the procedure to achieve economy. This 
concept is difficult to use for modem complex structuire for 
the same reasons which makes moment distribution obsolete. 

Phanjoo N. Ghista (js) also tried to design a 
full stressed structure, every section of which must reach 
its maximum allowable stress at some times during the appli- 
cation of alternate load systems. Obviously the cross- 
section is not kept constant. Initially the forces and 
moments in the structure are asstimed to satisfy the condi- 
tions of equilibrium between the external and internal forces 
and design parameters are obtained for full -stress condi- 
tioned. These do not satisfy the compatibility conditions 
and section properties are successively modified to satisfy 
compatibility. He could only solve two span beams and also 
recamnended the method for portal frames only. Hhanjoo N. 
Ghista (39) has also tried in another scheme to determine 
optimum frameworks under different load systems using the 
theorems of Mitchell and Maxwell. These methods along with 
mathematical programming methods axe only used for small 
stmotures. 

A 100 -storey building is designed by Pazlur R. 
Khan. {40) at ed. Here mahbers are proportioned by a IBM 
1620 and preliminary data for member sizes determined whicdx 
were fed to a bigger computer and the structure was analysed 



as a space structure. Louis i. Hill (41) designed a building 
frame using electrordo conputers. This was completely auto- 
matic beginning from initial sizing and dead load calculation 
to the final design. Elements were designed for optimum cost, 
taking into account the cost of concrete and forming. 

A language STRULL (42) (STEDctural Design Language) 
is developed for the design of structures. This language has 
the following operations: 

(1) Analysis Operation 

It can do preliminary and determinate analysis of 
a structure and then the stiffness and dynamic analysis of a 
structure. 

(2) Design Operation 

It can select member sizes and check these sections 
according to some codes of Specification. 

(5) Detailing Operation 

It can determine weight cost and can detail ccsm- 
ponents of a structure. 

Kinra and Penves (45) worked on a computer aided 
design. They analysed and checked concrete strsacturee. The type 
of structure studied by than is flat plate reinforced concrete 
buildings. The slab has a uniform thickness so also the gir- 
ders and beams which they designed according to ACI Code. The 
slab is idealised as girders (Tee beams) for the elastic 
analysis. All components are designed for both working stress 
and ultimate strength. 
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Beck (44) has developed a prograiDme to design 
columns of multi-storeyed "buildings subjected to bending about 
both the axis. In this programme he determines equilalent 
moments and for a given section satisfies the interaction 
formula. 

In this thesis the design of plane frame, is 
completely automated. The frame is analysed by the point 
iteration technique similar to that explained for trusses. 
Initial member parameters are determined by approximate ana- 
lysis. Contilever method is adopted for the analysis of Wind 
Load and for gravity load beams are assumed pined to the • 
columns to determine the column loads and completely fixed 
for beam moments. 

An \iniform I-section is designed for columns and 
beams. Since at the economic depth of I-seotions it has the 
minimum area> this area is achieved as fellows: Per a parti- 
cular depth the section is designed to satisfy the interaction 
relation. Then the depth is varied and the process repeated 
until the area of the section is a minimum. In short for 
each member an economic section, uniform throughout its span 
is designed to resist the forces it is subjected to. Then 
the structure is analysed for the above m^ber properties 
until the sections properties converge. Since the frame is 
made up of elements and elements designed are economical ones 
the frame is also an optimum one. 
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The thesis ocsritains the following chapters: 

Chapter II - Inalysis of Structures: 

The structures namely pin-jointed plane and space 
trusses and plane frames are analysed by point iteration of the 
stiffness equation. Computer programmes are developed which 
can analyse Beams, Plane Frames, Gables, Arches, Plane and 
Space Trusses. 

Chapter III - Feasibility of a Structural Design: 

In this section the fundamental equations of 
analysis are examined for both pin-jointed and rigid jointed 
structures. The feasibility of a full stress design for both 
cases are examined. 

Chapter IV - Desgin of Structures: 

In this chapter the preliminary analysis and 
structural design for both the type of structures are treated 
and suitable computer programmes are developed. 

Chapter V - Illustrations: 

Here seve3?al examples are illustrated. 

Chapter VI - Conclusions and Discussions: 

Conclusions and future work on the field is 

suggested. 



CHAPTER II 


STRESS AMD DEPOMATIOM ANALYSIS OP STEHCTUEES 

A structural analysis is concerned with the deteimi- 
nation of forces and deformations at every point of the stiruoture. 

In case of discrete structures like frames and trusses, interest 
is largely centred on the displacements of the joints and internal 
forces at the member ends. The jeint deformations and member end 
forces establishedthe ocanplete stress and deformation pattern in 
each member of a discrete structure. Once these forces and defor- 
mations have been found, the detailed calculation of the oonditionE 
at internal points of a member depends upon the characteristics i 
of the member. The material of the members are assiamed linearly 
elastic. For this type of structural members the internal forces 
and deformations are linear functions of the joint and members 
end actions. So for discrete structures the analysis may be i 

regarded as complete when the joint displacements and member end 

i 

forces are determined. The deformations and forces in a linear 
structure are interrelated and from one the other can be easily ; 

calculated. Accordingly there are two basic methods of analysis, i 

1 

The Pbrce method and the Displacement method. i 

In the Force method forces are taken as the unknowns 
and member deformations are calculated. The displacement compati- 
bility equations are satisfied frcam which unknown foTCes are 
evaluated. These are generally the redundant reactions. After 
this the stress or strain analysis of the entire structure is 
completed by the laws of statics and the stress strain relations. 



In the Displacement method 3oin.t defomations are 
taken as the tinknowns, member forces are calculated and joint 
equilibrium is satisfied from which joint deformations are 
determined. After this using statics and force defomaticn 
relationships the whole stress analysis of the structure is 
completed. 

Each method has its own merits. The Displacement 
method is more amenable to digital computers and hence this 
method is adopted. The final equations of the displacement 
method can be arrived at by several ways. Here slope deflection 
method is used to arrive at the equations of the Displacement 
method. 

2.1 Formulation of Slope Deflection Equations 

In this method member end deformations (the slopes 
and deflections) are the unknowns. The equations establish a 
relation which expresses the end forces of an elastic member as 
a f\inction of end deformations. This method may be used for the 
analysis of statically determinate or indeterminate frames or 
trusses. 

Here the equations are developed for a linearly 
elastic straight member having constant cross-section throughout 
its length. It is further assumed that the cross-section is 
doubly symmetric i.e. the centroid and the shear centre coincide 
resulting in no twisting moments for transverse loads acting 
through the oentroidal axis. 



Each meaiher is referred to an individual coordinate 
systan. In this S7’’stem the x-axis (also the axis-5) is taken 
along the centroid al-iaacie of the memher, y (also the axis 2) and 
z axis are along the principal axes of Inertia. This coordinate 
system is an orthogonal right-handed system and hence forth he 
referred to as the local system. 


2.2 Sign Ccnvention 

In the local system, forces and deformations are 
considered positive at a joint along positive coordinate axis. 


2.5 Plane FDMme Element 

For such an element there are three displacements 
along each memher ends namely the displacements, u, along x 
axis and v along y axis and the rotation, 0, about the z axis, 
and a compatible set of forces. 

y(2) 

^ x(l) 

Local Coordinate System 





Memher End Force 
u(l,l) 


u(l,5r*‘ 


I 


\ 


P(KI, 5) 


A ii(KI,2) 
^^^u(KI,l) 

^ u(KI,5) 


Memher End Deformation 


Pigt -f 



2*4 Nomenol ature 


A mem'ber is designated as I and KI, -where I refers 
to the left end and KI to the right end of the mamher. p(L,M), 
u(I,M) denote the force and defoimation at the end L and to-wards 
(or about in case of oouple and rotation) M axis respectively. 


The force and defomation vectors are respectively 
. P(I,2) 

! P(I,5) 


Moment positive about z-axis 
Shear force positive along y-axis 
Axial force positive along x-axis 


r 

u(l,2) 

u(l,5) 


Eotation positive about z-axis 

Sway defomation positive along y-axis 

Axial deformation positive along x-axis 


2*5 The slope deflection equations in the local system ares 
p(l,l) o 2E!I/L (2 u(I,1)+u(KI,1)) - 6EI/l(u(KI, 2) - u(l,2)) 
p(l,2) =. 6BI/L^ (u(I,i) +u(KI,i)) - 12EI/L^(4(KI,2) - ufc,2)) 
p(l,3) = -AE/L (u(KI,5)-u(I,5)) 

.. 2.1 

I, A, B and L refer respectively to the moment of inertia along 
the bending axis, the area of member, the Youngs modulus of the 
material and span of the member. 


Equation 2.1 in the matrix notation may be written 


as: 


'p(l,l)- W/l* 6 EI/l^ 0-' .u(I,T5| ^EI/C. -6EiA^ 0 :n(KI,l) 

i . ; ; 

p(l,2);- 'eEl/L^ 12EI/L^ O' Iu(I,2)+-6EI/L^^2BI/L^ 0 'u(KI,2)| 

0 0 ae/l: iu(i,3i . 0 0 - ae / l \ u(ki,3) 


_p(l,5^f 



In a condensed fom for fuirther work the above 
equation may he written ass 

p(l) = K(I,J) D(i) + K(J,I) ])(j) ,, 2.2 

where '^(1,1)1 !u(l,l)l ‘ iu(KI,ir 

! ' f 

p(l) ='p(l,2)i ]3(I) = !u(I,2) D(j) « u(KI,2) 

I . \ 

* I I 

p(l,5)i 'u(l»5V 

[Iei/l^ 6ei/l^ 0~' ' 2EI/L - 6 ei / l ^ 0 “1 

j j 

K(I,J) = 6EI/L^ 12EI/L^ 0 ; K(J,i) = i6ElA^ -12/EI/L^ 0 

t I 

j 0 0 ae/l ; 0 0 -ae/l', 

2,6 Truss Element 

For a pin jointed truss element may it he plane truss 
or space truss the members can carry only axial forces, hut each 
joint can have either two or three displacaaents along the coordi- 
nate axes depending on whether the tiruss is of the plane type or 
space type. There is no elastic bending •fnitoher. In the local 
system only axial displacement (i.e, along oentroidal axis) i.e. 
u(l,3) and u(KI,5) create forces in the member other displace- 
ments do not create member forces. Because of this the slope 
deflection equation for plane or space truss in the local system 
has the form: 

p(l,5) = - AE/l(u(KI, 3) - u(l,5)) 

p(l,2) = p(l,l) = 0 

2.5 

In the condensed form for future work it may be written as: 

p(l) = E(I,J) D(I) + K(J,I) B(J) .. 2.4 
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where p(l) = p(l,5); D(I) = u(l,5)*, D(J) » u(KI,5)j 

K(I,J) = AE/L; K(J,I) = -AE/L. 

Equation 2.4 is the slope deflection equation for a plane or 
space truss element in the Ir’cal systsa. 

2.7 Coordinate Transformation 

The above slope deflection equations are with 
reference to the local system. There are as many local ^sterns 
as there are members. In order to formiilate the equations of 
equilibrium all the above equations have to be referred to a 
common coordinate system. Such a system will henceforth be 
called as the Global System. This is also an orthogonal rectan- 
gular syston. 

In this section local vajciables will be transferred 
to global variables through a transformation matrix. 

Let X be a vector of local variables (may be the 
deformations or the forces in the Iccal system) and X be the 
same variables in the Global ^stem. Let T be the matrix that 
transforms the local variables to the global variables and be 
defined ass 

X = T X ..2.5 

T is known as the transfomation matrix. 

TUfhen one orthogonal set is transformed to another 
orthogonal set by rotation, the transformation is said to be 
orthogonal. Eor such a transformation (as is the case here) 



the* matrix T can he shown to have the property’s 

T””' s (inverse of ’T* matrix = Transpose of ‘T* matrix) 


2.8 Transformation Matrix for a Plane Prame Element 


y Y 



Similar Global vector be (6^ X y) 


Since Z axis does not change © is same in global or local 
systQa: 


X 


y 

e 

xy 


X CosiX + y SinO( 
X Sin(;(- Y CosX 

®xr 


In matrix notation: 

'"i 0 0 

0 -Cos A SinX 


© i 

xy! 


y 


! 

! 1 


1 0 Sino(Cosc< 


L.. 


©. 


XT 

y 

X 


or 


(t) 


I 1 
_ 0 
; 0 


0 c n 

I 

-CoscX SinX I 
SlnX CostX ' 


2.8 Transformation Matrix for Plane Truss Element 

It is same as plane frame element except that there 

is no © in this case. 
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I J ' 

x'l 

Y 


jCos^ Sino^! 
[Sinq-Coso^ i 

^Cos ^ Sino^ 
j Sino( -Cos ' 



Cos y{ 
Sin,3( 




; Cos:>^ 
! Sino< 


Sincx; , 

j 

-Cos.^ 

Sinc?( > 
-Coscx 


But since in local system we need only variables in 
direction only we shall only be using a part of the 
mat ion matrix namely; 


X = (Coso( Siu^) 


X ‘ and 


!Y 


ICosCK 


lY' 


' Kx) 

|SinO(i 


the X- 
transfor- 


Henceforth we shall define 
T = (CosO^ Sinfx) and T"^ = = 'cosf^i 

iSinCKi 


For the purpose of calculation T-matrix is always used and 
because of the above equations transpose of T is taken as its 
inverse although it has its origin with matrix. 


2.10 Transformation Matrix for a Space Truss Element 

It is similar to plane truss element except that 
the deformations or forces in Z direction (in global system) 
has to be oonsideared. 


form 


In this case the transformation matrix has the 


T = (ill) 
^ X y z'' 


« T 


ss 



where 1 , 1 , 1 axe the direction cosines of the local x~axis 
X y z 

in the global system. 


It may he noted that in the local system force 
deformation equations of plane frame or truss element has the 
form: 

p(l) = K(I,J) D(I) + K(J,I) D(j) .. 2,6 

where each variable refer to a particular matrix as the case 
may be. The transfomation matrix also transforms p(l), D(I), 
!D(J) from local to global system and has particular form for 
the particular case. The general case taking equation 2.6 and 
proper T-matrix will be dealt with. The particular oases mey 
be arrived at by proper substitutions. 


2.11 Equivalent Joint Loads 

A structure may be subjected to loads at the joints 
and loads along the member (only in case of plane frame elements). 
For forces acting in the span of the member, equilibrating joint 
forces (Fixed End reactions) are calculated. These loads when 
applied in the opposite sense form the equivalent joint loads 
in the local system. Local joint forces are transfonaed into 
the global forces and then added to the loads applied to the 
joints directly to form the total joint load. In case of truss 
members this case does not arises as there are no member forces. 
For a plane frame 


Let 


'S(I,1) 

jS(l,2)' 


be the Moment 
force in x and 

y direction respectively at a joint 



be the corresponding fixed end 
reactions due to member fomes. 


Let M(I,2) - 

I 

The total Joint Force is given by 

I r“ ” I ' 

P(l,2)l = 'S(I,2) j - t"^' + •;m(I,2)' 

F(i,5)J i,s(i,5)j 

.. 2.7 

2.12 J oint Equilibrium 

The internal fomes in the structure and the 
external loads must satisfy the conditions of equilibrium. 

As stated earlier for a discrete structure the equilibrium 
is satisfied when the joint equilibrium is established. For 
this the force-deformation relation is transferred to the 
global system and equilibrium equation at any joint is obtained 
by summing up the forces both internal and external along the 
axis of the global system, to zero. 

= 0 = 0 = 0 2.8 

d(l) = T D(l); d(j) = T D(J) p(l) = T P(l) .. 2.9 

Here variables in capitals refer to global system and variable 
in Knall case letters refer to local system. 

Substituting for local quantities in equation 2.6 
P(l) = T^ K(I,J) T D(I) + T^ K(J,I) T D(J). .. 2.10 

Substituting equation 2.10 in equation 2.8 

P(I) » (Z T^ K(I,J) t) D(I) + 2^*^ K(J,I) T D(j) 


2.11 



2.11 is the basic joint equilibrium equations. 

The summation A extends over all the members 
meeting at the Joint ‘I*; D(i) is the same for all members 
meeting at that joint where as D(J) differ from member to 
member. 

The equation 2.11 is a geneisl equation and holds 
good for both plane, space tirusses and plane frames In case 
of a plane truss it is a 2 x 2 matrix equation and for both 
spaae truss and plane frame it is a 5 x 3 matrix equation. 

One equation of the type 2.11 is written for each 
joint giving rise to the main system of equations, to be solved 
for the joint displacements. 

2.15 Boundary Equations 

It may so happen that for a boundary the joints 
may be restrained in one or more direction. A genei^ dis- 
cussion is given below for such joints. 

Let the Ith joint be restrained along 2nd direc- 
tion (l.e. in the direction of sway displacement). For this 
case: 

u(l,2) = 0; P(I,2) = 0 

Let r a> a^id fb'= ^t"^ K(J,I) T D 

Then the equation 2.11 may be written as: 

■P(I,1)^ fa,, a, 2 a,7 

,a,, ^2 a,3j |JI(I,5)J ';b(l,5); 
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Since P(l,2) = 0, this can be achieTsd by setting >= « 0 

^23 ^12 ” ^(l»2) = 0, can be achieved by setting 

®12 “ ®'22 “ ^32 “ is io say delete the proper rows and 

columns as shown and rewriting the matrix equation for the 


restrained joint we have: 


P(I,5): 


11 “13, 

I 

a,^ a- 


iu(i,i)! rbd.i)'^ 

+ 

,U(I,3), ;t>(l,3).' 


2.13 


~ 31 “33 _ _ 

Equations similar to 2.13 may be written for semi-restrained 
joints. If there are two restraints delete two rows and two 
columns. For joint which is restrained in all direction as 
can be seen there exist no equation for that joint * 


2.14 Evaluation of Member End Forces 

Once the joint displaconents ai*e obtained by 
solving the set of linear simultaneous eqviations the member 
end forces can be detemined by using the slope deflection 
equations. Since this equations are in the local coordinates 

f 

all global variables are to be transfoimed to local variables. 

The only variable needed is the displacements and can be found 
by relation 2.9« 

d(l) = T D(i); d(j) = T D(j) 

Substituting in equation 2.4 

p(l) = K(I,J) T D(I) + K(J,I) T D(j) ,, 2,14 

Equation 2.14 gives the member end forces in the local coordi- 
nate system. 



2.15 Method of Solution 


There are two methods for solving the linear 
simultaneous equations (2.1 1) for the joint displacanents. 

1 . Direct Method 

2. Iterative Method 


1. Direct Method 

In this case equations are written for all the 
joints. These equations in matirix notation have the fom: 


P = K D 

where P and D are respectively the joint force and 

joint deformation vectors. Both are 
(L X IT + m) X 1 - matrices 
K is the coefficient matrix (coefficient of 
joint displacements) also known as the stiff- 
ness matrix and it is a (L x N + M) x (L x N + M) 
matrix; where L is the total number of fully 
deformable joints and N denotes the number of 
possible joint displacements; L is equal to 5 
for both space truss and plane frame but is 2 
for a plane truss. 

M: Number of free displacements of all restraint 
joints (generally the boundary joints). 


2.1 6 Illustration 
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YlWr 01? THE STIPPITESS MRIX OF THE STRUCTURES SHORTS IN FIC-.2.2 
1 2 3 4 5 6 7 3 9 10 11 12 13 U 15 16 



: Five Band Stiffness I&trix 



Dimension of ’D’ matrix = 5 x (12) +(0+1+1 )2)= 40* 

Hence dimension of 'K' matrix is 40 3^ 40 • 

Joint equations: 1 throu^ 12. 

Here each joint has three equations. The equation for joint 1 
involves only the deformations of joint 1 along with that of 
2 and 5 only. The joints 2 and 5 will be referred to as the 
far ends and joint 1 itself as the near end joint. It can be 
easily seen that ary joint equation involve all far end joints 
deformations and the deformation of the joint in question. The 
maximum number of joint deformations among all the equations 
equals the number of bands in the 'K‘ matrix and the spread of 
all bands taken together is referred to as the band width. 

For the above structure by inspection the number 
of bands equal 5 sikL band width equals 5 x 9 “ 27» 

The figure F.3 gives the total stiffness matrix. 
Where a star (*) represent a 5 x 3 matrix oraas (x) 2x2 
matrix and dash (-) one by one matrix. 

Once 'K’ matrix is generated the dis|)lacements 
-1 

can be obtained as D = K F. 

Direct inversion is only possible for relatively 
p ma TT structures as the computer exhusts all memory for bigger 
ones. This is because lot maa^y null -matrices are to be stored. 
Hence this method is not suitable for large structures having 
many degrees of freedom. 



2. Iterative Method 


Advantage could be taken of the null matrices 
and storage can be considerably reduced by working with parti- 
tional matrices. In this case the total stiffness matrix 
is generated and stored and can be solved by several methods 
like guss Seidel, Jacobi iteration etc. In this case economy 
in storage may be achieved by storing only the bands provided 
it has a regular band pattern. It may be noted that for 
irregular structures and cl<^sed structures we do not have 
regular bands. In such cases storage economy cannot be easily 
achieved. For such cases the method developed below is helpful. 

Point Iterations 

In this method one goes from one joint to another. 
At every joint the particular joint equation 2.11 is generated 
and that joint deformation is calculated. The procedure is^ 
repeated for all joints several times until the iteration is 
stationary. This is similar to Gauss Seidel method# 

P(l) * K(I,J) T Ii(i) + K(J,l) T D(J) 

k(i,j) t d(i) = p(i) - k:(J,i) T 3)(J) 

Let A = K(I,J) T; B = K(J ,l) T 
A D(I) « P(I) - B I)(J) 

where both A, B, P(l) and D(j) are matrices. 

d(i) = a“^ (pCi) - B I)(J)) ..2.1 

Procedure: 

a) Number the joints in an optimal fashion. 
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b) Let JJ be the number of members meeting at a joint. 

c) For each member meeting at a joint generate 
T, K(I,J) and K(J,I) - matrices 

d) Per first cycle assume all joint displacements 
(may be assumed zero to start with). 

e) Calculate K(I,j) T and K(J,I) T D(j). 

Store them in matrices DENOM andUll’ER respectively. 

f) Detenaine the inverse of DBNOH 

g) Premultiply the Inverse of LMCM to (P(i)-TIPPER) 
thereby obtaining the displacements ])(l). Repeat 
the process for all joints. For second and sub- 
sequent cycles use the displacements of the pre- 
vious cycles and repeat for all joints several 
times until the joint displacements have stationaiy 
values. 

In fact this is the Gauss Seidel method. The 
oonvergenoy is quite fast. The oonvergency may still be 
accelerated by using overrelaxation factors. In this analysis 
the new displacements are obtained as: 

D(k) = I)(K) + 3(I)(k) - L(K-1)) .. 2.15 

where K refers to the iteration cycle and 5 the over-relaxation 

1 

factor. According to the above equation for the next iterative 
cycle the displacements of the previous cycle may be changed by 
a multiple of the difference between the displacements of two 


cycles. 



38 


The value of ^ generally lies between 0 and 
and a value of = (5.9 suits the framed structures. 

It is important to note that for the second design 
cycle the first cycles ^ deformations are taken and since the 
stiffness of the structure does not change radically only two 
to three more cycles for analysis may be enou^ to make the 
iteration stationary. 

The iteration is said to be stationary when the 
equilibrium of the structure is satisfied. Per this the total 
unbalanced load for each joint is detemined. This is added 
up for all the joints to achieve the total unbalanced loads 
for the structure. The iteration is said to be stationary 
when this total unbalanced load of the structure is negligibly 
small i.e. 2 to 5?^ of the external load. 
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CHAPTER III 

THE FEASIBILITY OP A STEHCTCEAL DESIGN 

The feasibility of the economic design of trusses 
and plane frames will be considered. The problan here is, given 
the geonetric configuration of a structure, subjected to a set 
of loads to detemine the design parameters to achieve an econo- 
mic design. 

Any valid design or analysis must satisfy the 
following fundamental equations. 

1 . Stress-strain Relation 

The structures designed are assumed to be in the 
elastic range. In the elastic range the stress-strain relation- 
ship is simple and involves a one to one correspondence of stress 
and strain. This means the stress and strain relation is independent 
dent of the loads the member has to cariy, i.e. the displacements 
and internal forces are linear functions of the external loads. 

Such a structure is called a linear structure. Since the changes 
in the geometry between the undeformed and deformed structure 
for a frame or a truss are really maall most of these structures 
under woricing loads behave approximately as linear structures. 

Because of this fact the stiffness and flexibility matrices 
are independent of the applied loads and the principles of 
superposition holdSgood. 
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5.2 Equilibrium Equations 

The forces (i.e. the reactions) at the ends of 
each member must be such as to keep the member in equilibrium. 
The reactions at the bo\indaries of the structure must be such 
as to keep the structure as a whole in equilibrium. The sum 
of the reactions acting at any joint must be equal to the 
external loads applied at the joint. This gives rise to the 
joint equilibrium, which are the main systsn of equations of 
the stiffness method of analysis as indicated in the earlier 
Chapter. 

5.3 The Conditions of Compatibility 

These refer to the displaconent compatibility 
equations although the above set of equilibrium equations may 
be called as the compatibility of forces. According to this 
the displacements of each of the member must be compatible with 
the displacsnents of the joints to which the member is attached 
or in other words the monbers so defoim that they can all fit 
together. These equations in the displacanent method is written 
as: 

T = Td •• 5*1 

where d is the number of independent joint displacements and 
V is the moaber end displacements of all the member of the 

i 

structure. T is a transformation matrix and only depends upon 
the geometry of the structure and also sometimes known as the 
geometric matrix. 
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For deleiminate structure the equations of equili- 
brium are enou^ for the ccaaplete stress analysis of the struc- 
ture and the forces in members are independent of relative 
stiffnesses of the members. The compatibility equations are 
mere identities in this case. The member areas can be designed 
according to any suitable procedure and may be designed to stress 
the individual mehers to their maximum capacity. 

For an indeterminate stiucture all the three cate- 
gories of equations are necessary for the stress analysis of the 
structure. The analysis and the feasibili'ty of a full-stress 
design procedure 110.11 be developed here for trusses and frames. 

3.4 Feasibility of a full -stress design for trusses 

A pin connected truss is a special type of structure 
and the load is carried by virtue of its axieO. stiffness. This 
stiffness is only a function of member areas and the length of the 
members. For a given goemetric configuration the variables are 
only the member areas and here the feasibility of a design is 
only with respect to areas of the moabers of a truss. As has 
al3ceajdy been pointed out the design of a determinate truss does 
not pose a problem because the areas of the members does not govem 
the analysis. But for an indeteminate structure the areas does 
influence the analysis. Hence the design of such structures are 
taken up here. 

For the truss, the elemental deflection is given 
by FEj/aE where A, L and E are respectively the areas, length 
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and elastic modulus of the material of the memher. The areas 
may be determined ass A « ^/cT » where the F is the member 
force andcris the allowable stress. F and cr remain constant 
throughout the length of the member. Deteimination of area 
by the above procedure satisfies only the strength crlijerion. 

It is generally uneconomical if stability criterion givejj,- 
a larger area of the mehers. Generally for usual columns 
the use of lattice and battened members give the same area* 
as determined by strength because the introduction of such 
members raise the buckling strength of the column to that 
value due to strength as in this case one has to look for 
local bucklings i.e. buckling of elements in between the 
lacings or battens and generally Eulers foimula for total 
buckling does not govern the design. 

For a full stress of in all members the elemen- 
tal deflections are given by 

■V « I cr= F/k •• 3*2 

£ 

Since C5", L and E are knoum for the inhole structure before any 
analysis or design, vs for all members are constants and known. 
According to the compatibility equations: 

d = •• 3»5 

In the above equation right hand side is all known as T is 
the transfomation matrix of the whole stiruoture and is a 
constant for any given configuration. This means before the 
analysis of the structure for full stress or equal stress in 
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all members, all the Joint deformations are constants and hence 
an an alysis is not at all neoessaryi 

This precisely means that for an indetsrmijiate 
truss because of compatibility restraints it is not at all 
possible to achieve full stress or same stress in all the 
members in the structure. In a determinate structure since 
the compatibility equations are only identities, such a con- 
clusion is not applicable. 

For more insight to these equations the compati- 
biliiy equations will be derived and then inteiJpreted using 
the principle of least woik or minimum strain energy or some- 
times also known as the Castigliano ’ s Second Theorems It may 
be stated as followss 

In any structuire the material of which is elastic 
and follows Hookes law and in which temperatuire is constant and 
supports unyielding, the first derivative of the Strain Energy 
with respect to any particular force is equal to the displace- 
ment at the point of application of that force in Ihe direction 
of its line of application. Mathanatically this theorem mey be 
written as: 

■^u/'c>Pn = Cn .. 5»4 

where U is the Strain Energy of the structure and Pn any load 
acting in the structure and ^n is the deflection at the point 
and along the direction of Pn. 
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5»5 Specialization of the Theorem for Solving 
Indeterminate Trusses 

Let designate the redundants, 'j* denotes the 
total nxamher of redundants and the final mamher forces. 

the member forces when unit forces are applied to the basic 
determinate structure in place of the redundant members, then 


U 11/2E (F^)^ 


3.5 


Summation is taken over all the members including the redundant 
members of the structure. 

T)U/r Rj == V "j •• 

where is the member forces due to external load in the 
basic determinate structure and are constants. 

.■^Fi/tRj » Ki ..3*8 

the summation drops out as, unless i «= J, "SF^/sR^ «» ^ 

K, is the member force when only one of the R.s act. 

1 J 

DU/-DR^ = Z(P. + I R^ L^A^E .. 5.9 

in the above expression the displacement along R- direction 

is included i. e. for the redundant members ’ j ’ P^ • and 

K . . >= K , =1 and we have / 

iD i ' 

which is the displacement along S^. There are equations 



of the type 'c)U/CR. and these are the oompatihility equations. 
The solution of 2)11/ = 0 gives the Redundant forces. 

where (T~^ is the stress in the members and 1*2 for tensile 
members and 1*1 for compression moabers. Substituting 
the fundamental equations becomes 

T cn (-1)^ K lVe. = 0} Let C. = L./E then 

•L 111 11 

cq = 0 .. 5.12 

where the sirmnation extends over all the members of the 
structure. 

Suppose all monbers are designed to achieve full 

stress then 0~ = fy * A constant and can be excluded from the 
i p 

svDBmation sign as tr- 0 and if B is also a constant 

P 

2^-0^ q q » 0 .. 5.13 

The above relation is a geometric relation and 
for any arbitary configuration need not be equal to zero or 
it is not possible to achieve full stress in all members of 
the structure. Since there are only 'j' such equations only 
'j* member stresses cannot be equal to the permissible values 
but are to be determined by the geometry. This makes the design 
of indeterminate trusses as complicated as its analysis. 

Illustration 

Suppose an indeterminate stinicture has m number 
of members and the redundancies be For a single load 



condition there axe '’j* equations in m unknown stresses. 

Because of this one can assign ar^ value to (m-j) member 
stresses and the ^3' stresses are to be obtained by solving 
the '3’ equations and these need not be equal to the penai- 
ssible stresses. Let the 'j* stresses obtained "by solving 
the above set of compatibility eqTiations be b^, in the ith 
member. If cr"^ is the permissible value then 
There are such equations to achieve full stress f^ = b^ 

or (-1 = 0 . But since this is the geometry of the 

structure, for achieving full stress the gemetry of the 
structure has to be predesigned. This is not the same as the 

. . 36 

‘Mitchell Stmctures ’ . This deals with the c«#figuration for a 
Statically Determinate Structure for a single load Condition, 
which has the minimum volume. 

All the above discussions are for one load condi- 
tion and the problem of two load condition is given below. 

For the same structure having m members and n 
redundants for two load conditions there are 2n compatibility 
equations althou;^ the number of redundant stresses are only 
n. In this case for the first load condition solve for the 
n member stresses such that (m-n) member stresses are equal 
to permissible value of the stress and the n member stresses 
are less than the permissible value. For the second load 
condition choose (m-2n) members and the 'n* member the stresses 
of which are not assxmied at the permissible value under first 



load condition and solve for the remaining 'n' member stresses 
such that these manber stresses remain within the safe limits. 
There are two difficulties encountered in the above procedure. 

1) The load conditions must be necessarily be two or 
more. For the same structure having m manbers and n redundants 
the design forces (higher of the forces under different load 
conditions) should be such that, at best, (m-n) manbers can 
have permissible stress under one load condition. Let 

^1i’ i 1 be the manber forces under the first load condition 

Pgi? i = 1 , m be the manber forces under second load condition 

i = 1 , m be the design forces in the members 

For all members fully-stressed at best (m-n) can belong 

to either or Pg^ separately. If say (m-n)+1 are from 

one of the load conditions then only n-1 compatibility equations 
are satisfied and one equation is not satisfied. Hence such a 
design will necessarily force one manber area equal to zero 
reducing the Indeterminancy of the structure to n-1. In short 
this means that in general even for multiple load conditions 
may not give a full stress design without forcing one or more 
areas to zero. For the multiple load conditions in which it 
is possible to design without violating the compatibility 
restraints will be referred to as the 'Real Multiple Load 
Condition* for futuire refeirence. 

2) Even thou^ apparently all the members of the 
structure are fully stressed under multiple load condition, 
a little thought will reveal that we are more or less 



designing a statically deteiminate structure for each load 
conditions separately. That is to say for each load condition 
a part of the stiruoture carry the load as a determinate struc- 
ture and little stresses are thrown to other manhers to satisfy?' 
the compatibility constraints. A different approach is shown 
for the design of such structure in the design Chapter. 

3.6 Feasibility of a Full-stress Design for Frames 
We start from the equation: 

V = Td 

In a plane frame the aim is to achieve full stress only at 
one section throughout its span. The member end deformations 
•v‘ depends upon the stress distribution throughout its span 
whidb are not equal to the permissible value. Hence forcing • 
one section to its full stress limit the displacements v are 
not determined and the compatibility equations are not 
violated. 

As was done before for trusses the compatibility 
equations from Castiglianos Second theorem may dbe derived as 

j«l 5=1 

N 

+ K (v. +T S. E.) S. ds/A.G. =0 * .. 5*14 

X —1^ X ^ X X X 

5=1 

The integration extends over all the members of the entire 
structure. are respectively the bending moment; 
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axial sforce and shear forces for the ith member when all 
redundants are removed, m^ and have the same meaning 
as and represent for unit moment and unit shear forces, 
ds is the length element A, E, G, I are the properties of 
the members and have the usual meaning and are defined earlier. 

The compatibility equations represent the defor- 
mation of a structure. For flexural members the deformations 
are predominantly due to flexure and direct stress or shear 
contributions to the deformations are ne^igible. For simpli- 
city only flexural deformation may be taken into account. This 
assumption reduces the above equations 3*14 io 
N 

+ Zj ®i ’“i “ 0 

N 

Let F. = M. + 27 where F. is the final moment at any 

X X i j X 

point in a structure. Then theabove equation reduces to 

F. m. ds/E^ = 0; F./l. = Cr/y. 

or Cr /y m. ds » 0 •. 5*15 

111 

In the above equation c5”^ is the stress at any point of a 
member. Our aim is only to fully stress only one section i.e. 
for one section of each member CJ". = restriction in 

no way effects the integral being not equal to zearo. The 
stresses at all other points of a monber can so adjustedastto 
make the integral equal to zero. 



50 


Even for one load conditions the compatihility 
constraints are not at all violated. This means for rigid 
jointed plane frames and also for space frames it is possible 
to full-stress one section of each member even tinder one load 
conditions. Although the analysis of such structures are 
complicated the design is not veiy different from that of 
a determinate structure except that there are member end 
forces . 
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CHAPTER lY 

DESIGN OF STEGCTUEES AND STEHCTDEAL CCSIPONENTS 

4*1 Approximate Analysis for Rigid Jointed Plane Frames 

For the exact analysis both the areas and moment 
of inertias of all the elanents of the stiructure are necessary. 
These parameters influence the convergency of the iterations. 

A better guess of these parameters may save considerable labour 
and computer time. Hence by some approximate method these para- 
meters are evaluated. 

A plane frame has to resist both lateral and gravity 
loads. Since it is a linear structure the influenee of lateral 
loads and gravity loads are determined by approximate methods 
separately and the solution superimposed to determine the total 
forces in the members. 

4.2 Analysis for Lateral Loads 

A tall structure is generally subjected to high 
wind loads. To simplify the analysis it will be assumed that 
the wind loads are concentrated at the floor levels. 

For this type of loads the structure may be assumed 
to behave as a Cantilever and the stresses are evaluated by 
ordinary beam theory. The method will be developed first and 
the assxmiptions will be made as and when necessary and finally 
they will be listed at the end of this section. 

For a Cantilevered beam the bending moment M and 
Shear force V for any section can be fomd out by statics. 



The stresses are evaluated by beam theory as: 

cr~ = Jiy/lj S » v/a .. 4*1 

where (y~±s the normal stress at a distance Y from the neutral 
axis, S is the average shear stress assxmied constant throu^out 
the section, A and I are respectively the area and the moment 
of inertia of the cross-section. 

The Cantilever analysis of a building flame is 
exactly similar to the Cantilevered beam. The moment ’M’ and 
the shear force V are the storey moment and the storey shear. 
Storey moment is the moment due to all loads acting above the 
storey level and storey shear is the algebraic sum of all loads 
above the storey level. Just as the moment varies from point 
to point along the Cantilevered beam so does the storey moments 
vary within a storey being maximum at storey bottom and minimum 
at the floor level of the storey. To get rid of this difficulty 
storey moments and storey shear forces are evaluated at the 
middle of the storey and assumed constant for whole of the 
storey. The normal stress and the shear stress of the Canti- 
levered beam becomes the axial stress and the shear stress of 
the columns of the storey under consideration. Since the 
shear stress distribution is assumed constant acixjss the 
cross-section of the Cantilevered beam the shear stress 
along the column section remain the same. If it is further 
assumed that the columns for a storey has some areas then 
the storey shear is equally divided among all the columns 
of the storey. 
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where and are respectively the shear force of a column 
at rth storey and number of columns of the storey. the 
column shear force remain constant for all the columns of 
the stor^. 

To determine the column axial force the moment 
of inertia of the column group of rth storey about their 
centroidal axis has to be evaluated. The centroidal axis 
passes through the centre of gravity of all the columns of 
the storey under consideration. If the colvimns have same 
areas as is assumed, then the centroidal axis passes throu^ 
the mid-point of the span of the storey under consideration 
if columns are symmetrically located about this point. 

= aE'R^ ..4.5 

where is the moment of inertia cf the rth storey and is 
the area of the columns of the same storey and assumed constant 
and equal to A. R^ is the distance of ith ooltunn of rth storey 
from the centre of gravity of the storey. 

= CT^ A = A MR^/A2;R. = *• 4*4 

where P^ is the cflumn axial force for ith column of rth storey. 

Since R^ is a constant for a storey the axial 
force varies as the distance famm the centroidal axis. This 
is analogous to nonaal stress in a beam which varies linearly 
from the neutral axis. 



The wind load is assumed to ax3t only at panel points 
and the span of the manbers are free frcan load. If the beam girders 
are assumed to be rigid compared to the columns, then the columns 
carries moments due to sway only. For such a case the column 
moments and beam moments change sign once along their span. Since 
the aim here is to determine the approximate values of the stress- 
resultants it may logically be assumed that the contraflexure 
points for wind load alone may lie at the centre of the beam span 
or column height. Once the hinge points are assxmied the stress 
analysis can be done by statics alone. 

CM = S X H/2; BM « j BS = •• 4*5 

where CM = The end moment of a col\mm 

S » Column s’^eax evaluated earlier 
H =« Height of column 
EM = End moment of a beam 

CM^ and CM^ ^ are respectively the column end moments 
attached to the beam from r and r-1 storey. 

BS = Shear force of a beam 

The summation sign is over all the columns of rtti storey. 

P, and P. - 6ire the column axial force for ith column of 
xr ir-1 

r atiH r-1 storey to which the beam belongs. 

Thus the stress analysis for the entire stMcture 
from wind loads is evaluated. 

4.3 Assumptions in the Cantilever Analysis 

The main assumptions may be listed thuss 



(a) Wind loads are concentrated at the floor levels 
only. 

(b) All colxJinns of a storey carry the storey shear 
according to their areas. If areas of all colinins 
are assmed equal the shear is equally divided 
among all the columns of the storey. 

(c) The axial force of a column varies directly as 
the distance from the centre of gravity of the 
colimin group of the storey. 

(d) There lies hinge points at the mid-height of 
columns and mid-span of each beam. 

4.4 Limitations 

This method may not be applied to a building 
frame where the beams are not continuous over all the columns 
of a storey i.e. column belongs to two stories. The method 
cannot be applied to building frames in which the column heists 
vary with in a storey as the storey moment and storey shear 
cannot be evaluated as before. But the method can be applied 
to any frame where the beams runs over all the columns having 
same height in a storey. There is no restriction for the number 
of columns in a storey. 

4.5 Gravity Load Analysis 

Since the beams are made up of mild steel and the 
section remains constant throughout the span» the sign of the 
bending moments axe not of importance to the designer except 
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for the possibility of some local instabilities. The magnitude 
of the design bending moment may lie somewhere in between 
simply supported bending moment and fixed end bending moment. 
For initial design the beams are designed for 120^ of fixed 
end moments. It is further assumed that the columns carry 
only the beam reactions and there is no distribution of fixed 
end moments in between beams and columns. 

The design load for any member is the addition 
of the absolute values of wind forces and gravity forces for 
the member. 

4.6 Determination of Approximate 'Values of the Displacemaits 
Although a free joint of a regular building frame 
has three degrees of Independent deformations, the axial de- 
formations are negligible when compared to the sway and 3?ota- 
tions except for the flower stories of a multi-storied build- 
ing. . The convergency of these displacements is not at all a 
problem and iterations may be stopped if the rotations and 
sways have stationary values. In this approximate analysis if 
axial deformations are assumed zero ’then the slope deflection 
equations for a member I-KI becomes? 

■P(I,1) i'^l/L 6EI/L^ ‘I jU(l,lT i2El/L -6EI/L^ I jTl(KI,l)[ 
'P(I,2)' j6El/L^12ElA^ ,n(l,2). ^EI/L^ -12EI/L^J _TJ(KI,2)J/ 

• ♦ 

It is proved in the numerical analysis textbooks 
that convergency of iterations of a set of eq.uations will 
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definitely occur if the absolute values of the coefficients 
on the main diagonal is greatei# than the sum of the absolute 
values of the other coefficients in a given row* 

It may be clearly noted that all structural 
analysis problems need not satisfy the above theorem# Hence 
although the problems will converge for stable structures 
the oonvergency may be extremely slow* To illustrate this 
two problems are considered below. 

Problem 1. (Page 180 and 176 - No. 64 and 63 respectively 
from Statically Indeterminate Structure by C*K* Wang) 


'a' 


(2) : B 

16 -(1) (i);0, 
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Pirobleni 1 

Since axial deformations are neglected the 
nxanber of independent displacements are: 

®B ®C 


as shoum in the diagram. The total stiffness matrix for 
the above order of displacement is: 


~6 

2 

1 

0 

-l/l6(-l) 

0 1 

0 ’ 

-2 

6 

0 

1 

-1/16(-1) 

0 . 


1 

0 

10 

2 

-1/16(~1) 

2 


0 

1 

2 

10 

-1/16(-1) 

4 . 

S , 

-1 

-1 

-1 

-1 

1/12(5/4) 

0 , 


' 0 

I 

0 -2 

where 

"4_. 

= a^/h 

0 

= Aj/i^ jzfg => 

5/12; 

A 2 ^2^. 

. 4.7 


Problem 2 


The deformation vector is: 




®0 ^ 

) 


shown in 

the diagram. 



: 46 
) 

15 

0 

-2/5 (-8) : 

• '! 

1 15 

cs 

80 

15 

-5/8 (-12.5)1 

, ®B ; 

1 

0 

15 

70 

-2.5 (-50) ; 

' : ®G i 

i -52 

-50 

-200 

19.9 (598) 

1 ^ (jz^/ 


^ = 


-d/20 



• • 


4.8 



A glance of the matrix A reveals that it has 
two extrsaely poor diagonal element namely a^^ The 

last two equations of ’A’ matrix in the classical analysis 
is known as the storey shear equations and these involve 
the sways deformations. 

The condition is not sc severe for the stiff- 
ness matrix of Problem 2. Here the diagonal elements 
are bigger than the off diagonal terns except for the shear 
equation (last equation in the 'B’ matrix). Although the 
convergence may be better in this case but still may be poor. 

Convergence for sway defoimations becomes 
extremely slow for single bay multi-storeyed buildings because 
of the shear equations. The convergency improves if the num- 
ber of bays increases but still may be poor. 

4.7 Improvement in the Convergence of Defoimations 

The convergence may be improved as follows* 

(a) Instead of if is assumed as the unknown as 
shown in the figure the coefficient in the *B* matrix of 
problCTi 2 becomes bigger and the diagonal terms fom a real 
strong diagonal improving convergency. But for problem 1 
the change in convergence may not be significant. 

(b) For a regular multi-storeyed building the sways 
or shear equations may be directly inverted and for arotations 
iteration may be resorted to. This may not load the computer 
memory as it may be assumed that each storey has only one sway. 
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(0) If a good guess may 6e done for the initial 
sways then also the convergence may he accelerated. These 
are taken from papers published in the early 20th century. 

Two of the methods are illustrated below. 

(1) Here the assumptions are: 

(1) Major portion of the sway is only due to wind 
load. 

(2) The floors are rigid. 

For such a case sway for each storey may be 

calculated as: 

M = 6 El or ^ = ML^/6EI •• 4*9 

Since moment for the columns are determined* and they are 
designed, can be found for any column of a storey and may 
be assumed constant for all the columns of the storey. The 
sway for the nth storey will be the sum of all the storey 
sways up to that storey. This may give exact results and are 
very useful for single bay multi-storey frames. Since the 
beams are assumed rigid the actual sways will become larger 
than the. above calculated sways 
209^ increase for all the sways. 

Illustration 10 ton - 


so it is better to assvane a 


■A 


B 


I 


/yy^ 

I 





Analysis by cantilever method gives Bending 
Moment, 3OO ton inch in columns. 

= ML^SEI = 0.5" if E = 15,000 
at 20^ increased sway becomes 0,6". 

The exact sway deformation even taking axial 
deformation « O.64". 

The stiffness matrix for the above problem for 
Q and becomes: 

.10 i/l -6 I/l^ ! 8.7 -0.036~' 

K = , P , or ’ 

i-6 i/l 12 ill? j 1-0.036 0.0005; 

where I/l = I25/144 = 0.8?. 

As can be seen it ihas very weak diagonal elanents 
and the sway displacement does not even converge after 50 cycles 
of iterations, but when initial sway as calculated above is intro- 
duced the convergence is achieved within four cycles. This 
method is adopted for the frame analysis because of its simpli- 
city. 

(ii) Here the assumptions are: 

(a) Axial defomations are! neglected. 

(b) Sway for a storey remains constant. 

(o) Joint rotations for a storey remains constant 
for all the joints of the storey, 

(d) To start with rotation for n-1 and n+1 storey 
is assumed same as the rotation of nth storey. 
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Procedure 

(1) Determine the external storey moment as was done 
in the Cantilever method. Let this moment be and for 
two consecutive storeys. 

( 2 ) Determine the internal moments of a storey. 

Manent for one column = M = = K(60 - 6 e) 

where 9 is rotation of the storey and R = A /h — sway angle 
of the storey. A is the sway and H is the height of the 
storey. 

Total moment of resistance of one complete 
storey = sxjmmation extends 

over all the columns of the storey. Similarly 
^ •• storey 

= Z' Kj, storey 

Adding and equating = 9^^ = 9 

= 6(ZK^ +2Kp 0 -6(7^1 +I \) •• A10 

( 5 ) Above is the storey equation. Now joint 
equilibrium equations become 

M = 0} summation extends over all members 
meeting at the joint. 


Equation for joint ‘A’ 

he 2®a * =®A - 5Sam ^AB - 5“b ^AB ^ ^ * 

®B * ^OA ®C ■ ° 
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This may be rewritten ass 

Wc ^ * V - 5(Ssm * "ab V * 


For the whole storey 

2®A(Z^=m l^Col’ -5' ®1 * V * 

^’'hA ®H ^ *M ®B ■'■ ^A ®0^ 


where T KL 

^ neam 

2.^Col 


K 


= twice the total stiffness of beams of a floor. 
= total stiffness of column of a storey 

V 

■= total stiffness of all columns of upper storey 
= same for lower storey 


% - ®N ■ ®o - ®p ■ ®q - ®r ■ ®1 
®B * ®B ° ®i' ■ ®a " ®j ■ ®1 ” ®T 

^2'^ ^ ^ ®B ’^A ®c> - ®]Z\ * ®r2^r * ®A/®Z^am 


Substituting the above simplifications in the joint equilibrium 
equations for all joints of the storey becomes 

“a * i^oi> * \Th * 

.. 4**12 


Now there are two equations for each storey namely 
^■Ih ®r ■*! “r V ®r 

h \ ('-SZ^Beam ^Z.^Ocl^ ®r * 

1 / 5(2 ®l^*Z®rV 
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This may be rewritten as: 

“a (■'aO * ■'aB ^ •'am) - * ®A3 ®AB^ 

% ■*■ ^ ®B ■'■ '^OA ®0^ “ ° 


For the whole storey 

®®a(Z"^^ I^col’ -5( Ih \ * I^r V * 

®M ■'■ ’‘m ®B ■'■ “^OA ®0^ 


where T K„ 

^ neam 

2^Col 


= twice the total stiffness of beams of a floor. 

= total stiffness of column of a storey 

\ 

= total stiffness of all columns of upper storey 
= same for lower storey 


®ll ■ ®N ■ ®o “ ®p ■ ®q ■ ®r ■ ®1 
®B - ®B ■ ®F ■ ®H - ®J ' ®L " ®r 

5,'=^ ®«*^®B*^®0> - ®lZ’^ * ®rZfr ®A/®Z\e=« 


Substituting the above simplifications in the joint equilibrium 
equations for all joints of the storey becomes 

“a * ®lZ^ * ®r^'^r ®1 

.. 4*12 

Now there are two equations for each storey namely 
^■Ih h \ -'h * ®r 

h \ - 2/5 (’-SZSeam *I‘'col> ®r * 

1/5(2 ®1 ^ ®r 
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Eliminating these two equations 

(“l + V/^ - ^ «r “ *l^olK * 

+ 1/3 (Z^i Sf 9j) 

.. 5.13 


If for the initial intention all 0 are assumed same and 

s 

equal to 0^ then 


8, - (M^+itp/eiKj^ 


earn 


4*14 


Let 2 K„.. 


Girder 


1/2ZS. 


earn 


then 


9„ - ..4.15 

Once 0Q is known can he found out from the storey moment 


equations. 


Pew iterations for each storey may given a fairly 
good solution. 

This method^ is due to John E. Goldberg (45 )• 




Fig. A ‘4 


— I'}' j'V- j 1 j^ j^ ~<*j 1^ 



lOTEEML MOIEHTS OF EESISm^CE OF A STOBE^ 
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4*8 Design of Plane Prame Members 

The ebject of a design is to provide members in 
a structure to satisfy the following requirementss 

(1) Seirviceability or functional requiremnets 

( 2 ) Economic requirements 

( 5 ) Strength and stability requirements 

4»9 Serviceability or Functional Eequirements 

The elements of a structure must safely discharge 
the function for which it may be built i.e. the serviceability 
criterion is satisfied. For example the water tank has to re- 
tain water i.e. the cracks in the tension zone of the concrete 
in the tank should not be of such magnitude which would allow 
leakage. Similarly the deformations of a beam supporting a 
glass panel must be small enou^ to safely support over head 
panel. 

The serviceability of a plaue frsme generally 
refers to the defoimation which should be within a permissible 
value satisfying the functional requirements. The permissible 
deflection of a beam is usually relative to the span but some- 
times it may be independent of the span if it has to discharge 
some special requirements. For the above exampls of the glass 
panel over the beam the deflection has to be such that no 
damage is caused to the structure above and hence it has to be 
independent of the span. But these are special cases and the 
design here is governed by the general case. 



The Indian Standards Specification (IS-SOO) limits 
the maximum deflection to span/325. 

The correct determination oi the deflection is 
as difficiilt as the analysis itself. Instead deflection will 
be related to the depth of a simply suppo3rted beam with an 
uniformly distributed load throughout its span. Since the 
deflection of the members a rigid frame is expected to be less 
than the simply supported beam this calculation keeps the 
design on the safe side. 

(5 IiL^)/(48 El); » Md/21 

(51L^)/24 Ed 

» 1/325? = 20 ksi; E = 30,000 ksi 

L/d = 22.2 ..4*13 

where M « wL /s? w being the uniformly distributed load over 

the span L; ^ is the central deflection; and d is the 
max 

depth of the beam; E,I has the usual meaning. 

If the depth of the beam is at beast l/22«2 of 
its span the criterion of serviceability is satisfied for a 
mild steel beam. Since the economic depth of the beam gene- 
rally lies in between l/l5th to l/lOth of its span service- 
ability is not at all a problan fer mild steel beams. 

4.10 Economic Requirements 

According to this requirement, for a given bending 
mcmient, shear force, axial force and span, the section designed 


-^max 


max 


9 

? mas'^^ 



is economical when the area of the section is a minimim. 

Instead of volume, area is taken as the criterion because 

« 

the section is kept, constant throughout the span of the 
member. An I-section is designed for the elements of a plane 
frame. For such an element initially the area of the section 
reduces as the depth of the section increases from service- 
ability depth to a certain depth and thereafter because of 
local and total instabilities of the member, the area starts 
increasing for any increase of depth. The depth at which the 
area is a minimum is temed as the economical depth. The aim 
here is to achieve this area. 

For a particular depth hi^er than the service- 
ability depth the area of the I-section is deteimined which 
satisfies the strength and stability criteria. The above pro- 
cedure is repeated for several depths until the economic depth 
is attained* This gives the minimum area #f the member. 

4*11 Strength and Stability Requirements 

This requirement is, in fact, the prime criterion 
for the design of a structural element of a building frame. In 
building frames, the structural members are generally classified 
as beams and columns according to their position and type of 
the predominant loads. In the general frame design such a 
classification is not necessary, instead a general formulation 
is made to design an uniform I-section throughout the span to 
resist the given forces and automatically either a beam or a 
column is designed according to the forces the section has to 


carry. 



69 


For a particular depth, d, a uniform I-section 
is designed to satisfy the strength and s.tability criteria. 

The problem here is to design an I-section for 
a given set of loads namely the bending momoat, the axial 
force and the shear force acting on a span of length L. Each 
elements of the I-section namely the Web and the Flanges are 
designed for strength and stability. 

4.12 Strength Criteria 

For an element in direct compression and bending 
about one axis the criteria of strength is governed by the 
interaction formula: 

f /F + f,/F, C 1 .. 4.14 

a' a b b ^ 

where f and f^ axe respectively the axial and bending stresses 

Of U 

of the section and F^, F^ are respectively the peimissible 
values for the same “Stresses. The interaction formula is satis- 
fied by iterrating on the area of the section. 

For the initial determination of f and f, it is 

Qi D 

assumed that the flanges carry the total bending moment and 

2/5rd of the axial forces, f and f, are, proportional to the 

bending and axial forces in the flange; and 

F = F, = , 

a 0 

being the permissible value based on yield criterion. 

Since F = F, the iteraction formula becomes 
a D 





1 '^ 

Force in each flange due to moment = M/d 

Force in each flange due to axial force = P/5 

(M/d)/(P/3) = K * f^/f^5 f^ = 5M d 

^a “ ^a d) = CT- 

- P, = (C" P d)/(P d + 3 M) .. 4.15 

* ( 3 M<r)/(P d + 3M) .. 4.16 

The above formulae are necessary for a preliminaiy value of 

the section area. The initial area is higher of the two values 

A = P/f or A = (1.2 M/d + P)/ 

3i 

where F is the direct compressive force. 

Since the section designed is an I-section the 
approximate moment of inertia is given by 

I = 0.38 A X d^; Z = 2l/d 

f^ = P/A? f^ = m/Z 

Let BETA = f, /P, + f /P 

b' b a' a 

The value of BETA determines the safety of the 
section against strength. The structure is safe or unsafe 
according to whether BETA is higher or lower than one. In 
both cases if area is increased BETA times the iteraction 
fomula is identically equals to unity and the section is 
safe. 

4.13 Design of the Web Plate 

The section designed is an unifom I-section 
without any vertical or horizontal stiffness. The web 



71 


of the I-section has to he checked for buckling due to shear, 
bending and axial compression. 

The critical stress at which buckling is imminent 
for a rectangular plate with various boundary conditions and 
various load conditions for load in its own plane is given by 
the formula 

S = (K E)/(l2(l-'J^)(b/t)^) * .. 4.17 

where K = a constant dependent upon boundary conditions and 
load conditions (shear, bending or axial load). E and are 
Youngs modulus and Pinssons ratio; t = thickness of the 
plate; b = length of the loaded edge of the plate, except 
when the plate is subjected to pure shear it is the 
smallest lateral dimension. 

Substituting for E = 50*000; = 0.3 

= K 27100/(b/t)^ .. 4.18 

4.14 Determination of the Thickness of the Web Plate 
To Prevent Shear Buckling 

In this C5ise the boundary is taken to be simply 

suppoarted along all edges and a/b or span/depth ^* 01 “ 

the above conditions the value of K is a minimum. 

*This formula along with other formulae of this Chapter 
are quoted from the design of steel structures by 
Gaylord and Gaylord. 



For this case the K value for pure shear load is K = 5»3 
upon substitution: 


= 144000/(d/t)‘ 

cx 


4.19 


Allowing a factory safety of 1.66 against buckling 


al 


87000/(d/t)‘ 


Hence the allowable stress is either permissible 

value of shear stress based on yield criterion for mild steel. 
For safe design the actual shear stress should at best be 
equal to the allowable shear stress. 

1.25 v/dt^ 


max 


= (1.25 X V X d/87000) 


1/5 


.. 4.20 


The factor 1.25 is used because the shear stress v/dt^is the 
average value, whereas the maximum shear stress for an I- 
section is generally 1.25 times the average shear stress. If 
the web has a thickness t^ it can carry safely the shear force 
V without buckling in shear. 


4.15 Buckling due to Bending and Axial Force 

The section is subjected to normal stress due 
to both bending moment and axial compression. This normal 
stress due to both the effects must lie within the permissible 
limits. 


The web plate of the I-seotion can buckle either 
in bending or in direct compression. The thicknesses for the 



web plate are detennined for both the cases separately so as 
to achieve permissible stress based on yield criterion with- 
out buckling. Since the pemissible value for mild steel is 
same both in bending and axial stress, the design thickness 
is taken as the sm of the two thicknesses. 

i) Bending Buckling 

Substituting suitable value for K the allowable 
stress in bending for a web plate to prevent buckling is given 

®cr “ 650,000/(h/t)^ *. 4*21 

If a factor of safety of 1.53 is taken against buckling due 
to bending and maximum allowable stress is taken as 20 ksi 

(h/t) » 175 “ 1^/173 4*22 

Hence if the depth is within 173 times the thickness of the 
web bending buckling does not occur and maximum allowable 
stress of 20 ksi may be achieved with a factor of safety of 
1 . 55 . A low factor of safety of 1*55 is assumed because the 
< K» value is taken for a plate with all sides free^ and this 
is a conservative value for the web of an I-section# 

ii) Compression Buckling 

Substituting proper value of the allowable 
stress i*e. comp 3 ?ession to prevent buckling due to compression 
beobmess 

S « 204000/(h/t)^ •• 4*25 



■ » 
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To achieve a 20 ksi stress in compression with a factor of 
safely of 1.55 

h/t = 87 or t = h/87 •• 4*24 

c 

t, and t are the two hounds for the weh thicknesses for pure 
D c 

compression cr pure handing. In the actual case since hoth 
bending and axial forces act, the thickness for hoth the cases 
may he determined substituting the actual stress in the weh 
due to bending and due to compression. 

h/t^ . or - h(f^)''/^/390 

and . h(fj,)'’/V775 •• 4.25 

where f and f, are respectively the actual stresses in the 
CD 

weh which are calculated previously. 

The thickness to prevent buckling due to bending 
or due to axial compression is avoided by taking the sum total 
of the two thicknesses as the allowable value for mild steel 
section due to direct compression or bending 'is the same. 

The above foimulae are derived to prevent buck- 
ling due to normal stresses or shear stresses. 

In a rigid jointed frame the marbers are gaierally 
subjected to hi^ moments and concentrated shear forces at the 
ends. These give rise to combined shear and direct stresses. 
For such cases it is extremely difficult to specify the exact 
thickness of the web to avoid buckling due to the ocSibined 
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action of the shear stress and normal stresses. Since the thick- 
ness to prevent shear buckling is usually small and is gene- 
rally less than half the minimum thickness required by the I.S. 
Code of Specifications (which specifies 0.25" as the minimum) 
t^ will be taken as 2t^ . In this case the web shear stress 
will become 

v/v^^^ = 1/0.25 or V = v^^^ X 0.25 

The interaction formula when both normal and 
shear stresses act is given by 

(v/vcj,t)^ + ..44*26 

where s and v are the bending and shear stress of the section, 

s . , V . , the allowable values for the same. Substituting 
crt ’ crt ’ 

■^/^crt “ s/s^^ = 0.95. 

providing twice the thickness for shear re- 
quirements, the normal stress only has to be lowered by 55^ 
gnri for all practical purpose if the ihickness due to normal 
stress is determined assuming the shear force zero the section 
designed is safe. 

4.16 Design of Flanges 

The flanges carry a major part of the bending 
moment and axial force. The area of the flange may approxi- 
mately be determined as 

= (l/<r ) (M/d + f/ 5) - 1/3 \ 


.. 4.27 
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In the above expression the bracketed term refers to the 
approximate force the flanges may have to carry if the web 
does not carry ary bending moment and l/5 Axial force 
•F’ is being carried by the flanges. But since web carry a 
part of the bending moment, one third the area of the web is 
deducted as the contribution of the web for carrying the bend- 
ing moment. In the textbooks the web contribution to resist 
bending moment is generally taken as one sixth the web area 
but for shallow depths and for low bending moments the web 
contribution may be as hi^ as I/3 itsarea. 

4.17 Buckling of Flanges 

Both direct load and bending moment causes direct 
stresses in the flanges of an I-section. The compression 
flange may buckle due to the combined conpressloa. The. criti- 
cal compressive strength of a bean flange is given by 

CTc “ (nd)(E ly GJ)^/^)/(2I^ L) ..4.28 

where I^ = 2 btV5> 3!^ = tbV6; J = 2btV5 

Substituting for E *= ksi and G = 12000 ksi for mild 

steel and with a factor of safety of 1.66. 

CT^ = 12,000/(Ld/bt) .. 4.29 

where L, d are the span and depth of the beamj b, t are the 
breadth and thickness of the compression flange. 

In determing the area of the flange by equation 
4.27 proper value for is substituted. This value is decided 
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on the basis of yield or on the basis of buckling due to 
compressive stresses. The lower value of is used. This 
keeps the section in the safe side. 

Once the flange and web elements are designed 
the exact moment of inertia I, the area A and the radius of 
gyration r are calculated. 

4.18 Compression Buckling of the Section as a Whole 

The critical stress in a column is generally 
eacpressed by the Euler's formula as 

^rt ..4.50 

The parameter is substituted to take care of 

the proper boundaries. For non-sway cases the value of L is 

0 

generally smaller than the clear span of ■ftie column. But for 
sway case it depends upon the factor 

G = (iyL^)/(lg/Lg) (46) ,. 4.31 

where the subscript c refers to the column and g to the beams. 

In a frame the girders are heavy and since a part of the slab 
over the girder usually act in combination the stiffness value 
of it is extremely high. Due to this fact the factor G beccmes 
extremely small for usual building frames. For G equal to infinity 
the effective length becomes twice "the clear length and for 
G equal to zero L becomes the clear span. In our case since 
G is small is taken as the clear span. This is also as 
per the I.S.I. 800 Specifications. 
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Finally the interaction formula is satisfied 

(P/l)/^^^ + (M/zVc-^ 

must be equal "to 1 within 5?^ accuracy. If the interaction 
formula is not satisfied within 59^ tolerance, the area is 
either increased or decreased in steps of few percents of 
the total area and the whole process repeated until the 
relation is satisfied. 

This gives a section for a particular depth. As 
already mentioned under economcc requirements the process is 
repeated for higher depths until the minimum area of the 
section is achieved. 

4.19 The Design of a Section may be summarised as follows* 

i) Choose a depth to satisfy the serviceability 

requirements. Let this depth be ’d' = span/22, 

ii) Determine the thickness of the web according 
to the equations. 

4.20, 4.22, 4.24, 4*25 and satisfying I.S. 
Specifications for minimum conditions, 

iii) Determine the area of the flange from equation 

4.27. Assume a thickness (minimum 0,25) and find 
the spread of the flange. This spread must not 
be greater than I6 times the thickness (Code 
requirements I.S. 8OO) 

iv) Determine actual area and moment of inertia 


of the section. 
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v) If ’BETA* is hi^er than 1 the section is unsafe 
and increase the area of the section. Distribute 
this increase in area in the web or flanges accor- 
ding to their areas. 

If BETA is less than 1 the section is oversafe. 

In this case the process is repeated and checked 
if any of the areas of the flange or the web can 
be reduced without violating any of the buckling 
or strength constraints. 

vi) Repeat the whole procedure for higher depths 
until the area of the section has an absolute 


minimum value 
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4.20 Design of Trusses 

A preliminary design is necessary to accelerate 
the rate of convergence in the analysis and in the design. 

The convergency of the analysis depends upon the magnitude 
of the diagonal elements of the stiffness matrix. When all 
the areas of the memhers of a truss are p.ssumed equal, the 
relative magnitudes of the diagonal elements of the stiffness 
matrix remains unaltered although the total stiffness of the 
structure may change which in turn may affect the deflections 
of the structure. In general, in the stiffness matrix of a 
truss, the diagonal elements are the predominant elements in 
any row. This is because the manber stiffness matrix of a 
truss has squares of the direction cosines as the diagonal 
elements which are always positive and the equilibrium equa- 
tions of a joint which lead tc the structure stiffness matrix, 
adds up all the positive terms for all the maibers meeting 
at the joint. These elements appear on the diagonal of the 
stiffness matrix giving rise to a very strong diagonal. 
Generally the off diagonal terms are single stiffness elements 
and are smaller . This clearly shows that for trusses the 
choice of member areas dc not influence the convergence for 
analysis. The trusses generally have low degress of indeter- 
minacy as the addition of one member only increases the 
degree of redundancy by one unlike in frames where the 
addition of a member increases the indeterminacy by three. 
Since changes in the member forces depend mostly upon the 
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degrees of indeterminacy .• for a truss the member areas 
do not influence the m^ber forces to any appreciable extent 
like the frames. Because of this any arbitaiy areas may be 
assumed for the members as preliminary areas. The volume of 
the structure obtained in the first cycle of design and the 
volume when all members are fully-stressed does not differ 
much. 

Because of the above facts most of the trusses 
do not need a preliminary design. In the ccaaputer programmes 
developed equal areas for all members are assumed to start 
with. 

If a preliminary design becomes necessary, then 
a statically determinate truss may be chosen out of the in- 
dete'minate structure. Member areas m^ be determined for 
the determinate structure and the volume of the truss may be 
calculated. On the assumption that an indeterminate truss is 
lighter than the determinate one the above volume or few per- 
cent less than the volume may be distributed among all the 
members of the indeterminate truss according to the reciprocal 
of monber lengths , since the forces carried by a member 
depends upon the reciprocal of its length. 

For the final design the trusses may be classi- 


fied 'into three categories. 
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4.21 Detenninate Trusses 

The equations of statics aie enou^ for the ana- 
lysis of these structures. The member forses axe independent 
of the member areas and the preliminary design is the final 
design. The members may be proportioned to achieve full— stress 
condition either under one or multiple load conditions. 

4.22 Internally Indeterminate Trusses 

Fo"^ a structure of this type the external reaction?* 
can be determined by the Laws of Statics alone. Since the stiuc-* 
ture is a determinate one for the external reactions but the 
internal member forces are indeterminable by statics these 
type of trusses may be called as internally indeterminate 
trusses. 

■“'or a plane truss three independent equations of 
statics are available to determine the extemal r-?actions and 
if the reaction components are three then it comes under this 
category. For a space truss these reaction ccanp clients are six. 
The above number of reaction components make each truss stable 
and externally determinate. 

Fig. 5 shows an internally indeterminate truss. 

Although it is an indeterminate structure the 
member forces, for the determinate panels namely panels 
1 - 2 - 5 - 4, 5 - 7 - 6 - 8, 8 - 10 - 11 - 12, and 

■11 -12-15-14 can be determined by statics. These forces 
remain independent of any member areas of the structure. 





For example! taking a section in the panel 
5 “ 7 " 6 - 8 the member forces can be detemined by the method 
of section as follows! 

Taking moment about the joint '8' 

^57 “ - ^ ^ ••4.33 

Taking moment about jf int ' 5 ’ 

^68 “ ^12 (a + b)/d .. 4.54 

Taking moment centre at joint *6* 

F 53 = - + b) + X d)/d .. 4.35 

Similaa^f oroes in the other determinate panels may be evaluated 
by the equations of statics. 

An internally indeterminate truss may be divided 
into two Complexes. 

(a) The Determinate Complex 

This includes all the determinate panels. The 
member forces of this complex only depend upon the external 
loading and are determinable by statics alone like any deter- 
minate structure. The member areas of this ccanplex does not 
influence the force of any monbers of the structure. These 
members may be designed just like statically determinate 
trusses. The member proportions should be kept constant 
throughout the design of the structure. This may save lot 
of human and computer time than designing all the members 
of the truss for each design cycle. 

(b) The Indeterminate Complex 

This includes all the indeterminate panels. 



To evaluate the foioes of members b^onging to this complex 
indeterminate analysis has to be canded out. The forces 
of members belonging to this complex are influenced by areas 
of this complex alone. This can be readily arrived at by 
strain energy approach. 

Let TJ = + Ug .. 4.56 

where U = Total Strain Energy of the whole stiructures 

, Ug, the strain energy of the determinate and 
indeterminate Ofmplexes respectively. 

The fundamental equations of Force method ares 
dU/ dR^ S3 r^ ; where is the iredundant force and 3^ is the 

displacement along R^ 

But dU/dR^ = dUg/dR^ as dtJ^/dR^ = 0 as TJg is ca constant and 

is independent #f R^ as member forces of this complex axe 
statically determinable. 

dU^/dR. = r. ..4.57 

or the strain energy of the indeterminate complex alone deter- 
mines the redundants. 

This complex can separately be dealt and me®ber 
forces determined. Since the member forces are a function of 
the member areas of this complex alone at each stage of the 
design these areas are to be altered till full-stress in all 
members or minimum volume is achieved. 

Many of the bridge trusses contain a few indeter- 
minate panels and cone imder this category. If the above pro- 
cedure is adopted designer can complete the design in »uch 



less a time than necessary for the usual procedure in which 
no daaarcation is made about the complexes. 

4.25 Externally Indeterminate Trusses 

A truss is defined to be externally indeterminate 
whCT. the external reaction components cannot be determined by 
the laws of statics alone. If the reaction component exceeds 
three for a plane tiruss or six for a space truss then it comes 
under this category. These reaction components are not only 
dependent on the external load but also are influenced by the 
stiffness of all the members of the structure. Since forces 
on all members depend upon the external reactions member 
forces are influenced by the stiffness of all the members. 

The proportioning of all the members has to be carried out 
at each design cycle. 

4.24 Method of Design 

In this procedure only indeterminate Complex 
members in case of an internally indeterminate truss and all 
members of an externally indeterminate truss are designed. 

Other members may be proportioned to achieve full-stress 3 ^st 
like a determinate truss. 

As earlier indicated there are two separate types 
for the design of indeterminate trusses. 

1 . Design for Multiple Load Condition 

2 . Design for Single Load Condition 

4.25 Design for Multiple Load Condition 

It is earlier pointed out that even for multiple 
load condition, the conpatibility equations may not allow all 



the mauhers to be fully stressed. Althou^ this is true, this 
becomes a rare case. This is because Ihe real multiple load 
condition (defined in the Chapter - feasibility of a full-stress 
design) depends upon the areas of the members of a truss and 
these areas keep on changing at each design cycle. So even 
though the area distribution initially may not give a real 
multiple load condition, the areas can distribute in such a 
manner that at a subseq,uent design cycle real multiple load 
condition may be achieved. Generally except for some "patho- 
logical” cases full stress in all members may be achieved for 
a multiple load condition. 

The design for a truss under the above irestric- 
tions may be broken up to the following steps. 

(1) For an assumed area of the members, the member 
forces may be evaluated. 

(2) For the above analysis proportion the areas 
of the membera to achieve full stress in all 
the members. In this design, for compression 
members, compression buckling has to be pre- 
vented. The design SUBEODTINE developed for 
the design of plane frame maubers may be used 
if moment and shear forces are equated to zero. 
The use of lattrice and battened manbers make 
it possible to avoid the compression buckling 
altogether and permissible stresses may be 
achieved in most of the cases. 



(j) For the above areas the analysis may be carried 
out onee again and member stresses evaluated. 

The absolute maximum of these stresses is cal- 
culated and determine the factor as* 

The absolute mflYlmuTn stress in any monbers of the stiruoture 



The permissible stress 

.. 4*58 

If 1 then scsme of the members are overstressed 
and for p<1 all member stresses are within the 
allowable value and it is a feasible design. 

( 4 ) For p 71 prori'te all the areas of the members by 
the factor jg . By this process the analysis re- 
mains unchanged since the relative stiffnesses 

of the members remain the same. This isnow a 
feasible design. 

( 5 ) For redesign the members to achieve full 
stress. 

The design is said to be complete when the volume 
for the feasible designs of steps 4 and 5 has a minimum value 
or remains constant as the design iteration proceeds. 

If the volume for the feasible designs are plotted 
against the design iteration it is found that either the volume 
reduces till it remains a constant or it reduces and then in- 
creases (see Illustrative problaas) and finally remains a cons- 
tant as iteration proceeds. In the case where the volume 
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reduction is monotonic the full-stress design is also the 
min i mum weight design otherwise the full stress design is not 
the minimum wei^t design. 

For the above pTOcedure no mathematical proof 
is foTind as yet but the plot of volume vs iteration for some 
problems gives lesser volume than that of a full-stress design. 
It should not be thought that full-stress design has approxi- 
mately the minimum volume as even for a lesser volume than 
this the relative member areas may be so distributed that 
the member forces and stresses may be lower than the fully 
stress values. This may be the real absolute volume of the 
structure. 

4.26 Design for Single Load Condition 

46 

In 1900 Celley pointed out that for one load 
condition generally a deteminate structure has the minimum 
weight. For this wei^t the deflections are also a minimum. 

It was also pointed out in the Chapter on the Feasibility of 
a fully-stressed design that a fully-stressed statically in- 
determinate structvire can not exist for one load condition 
as it violates the equations of compatibility. So it is not 
at all desirable to go for an indeterminate truss to support 
only one set of loads. .But if it is necessary to design a 
statically indeterminate stiucture may be to minimise the 
deflections then the same procedure as for multiple load 
condition may be adopted except that instead of proportioning 
all the members to attain full-stress limits only M-N members 
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be designed to achieve full stress, where M is the total number 
of members and N is the number of indeteiminance of the structure. 

By doing so the equations of compatibility are not violated. 

47 

Wilber and Norrice have pointed out that for one load condition 
if a statically indeterminate structure be designed to fully 
stress all the members at each cycle of design then finally 
several member areas will reduce to zero to yield a statically 
determinate structures. This is because compatibility laws 
are violated. But the recommended design does not give a 
determinate structure as all the laws of analysis are satisfied. 

The choice of (M-N) members out of the M manbers is 
is a really difficult problem. This choice may decide the volume 
of the structure. This (M-N) quantities may be chosen according 
to the member stresses or member forces. Suppose an indeter- 
minate truss has 15 members and the redvindants are five then 
either choose the 10 highest member forces and proportion 
these members to achieve full stress or choose 10 highest 
member stresses and design these members for full-stress. The 
volume of the structure nay be determined for both the cases 
and the lighter structure may be chosen for economy. This 
volume will definitely be more than the volume for a deter- 
minate structure as essentially we are designing a statically 
determinate structure out of the total structure and other 
manbers of the structure carry little load only due to the 
laws of oompatibilily. This structure may have lesser def- 
lections than the determinate one. 
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In conclusion, it may be pointed out that we are 
intiroducing arbitrary means to satisiy the laws of analysis 
and there is no constraints to minimize the volume at all. So 
this design may not have the absolute minimum volume? 



ILLUSTRATIVE PROBLEM NO 
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CHAPTER ¥I 

STBfflim OP RESULTS, CONCLUSIONS iND DISCUSSIONS 

Since the process of resisting external loads 
for rigid jointed and pin- jointed structures are different, 
these structures are separately treated to hi^xLighten the 
important features involved in their design. 

6.1 Pin- jointed structures 

Reference t Chapter IT : Illustrative Prohlon No. 1 

This is a plane truss with one degree of indeter- 
minacy and is subjected to two load conditions. It took six 
cycles to achieve the stationary values of the joint displace- 
ments by the point iteration teohni(lue wi1ii an initial 
assumed zero displacements for all joints, ill the member 
areas were assumed to be equal at the outset. The first 
feasible solution was obtained by prorating all the areas of 
the members in aocordanoe with the maximimi member stress. 

This volume was 567.48 cubic inches. The members were then 
proportioned to achieve full stress for all members in either 
of the load conditions. For these member areas and joint dis- 
placements, only three tycles of iterations were required for 
the point iteration procedure to converge. The member stresses 
detemined and member areas prorated to achieve another feas- 
ible solution. Graph 5*1 shows the volume for feasible struc- 
tures versus the cycles of iteration. The minimum volume 

sixth iteration. For this volume only 


was obtained at the 
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CHAPTER YI 

SUMMARY OP RESULTS, CONCLUSIONS AND DISCUSSIONS 

Since the process of resisting external loads 
for rigid jointed and pin- jointed structures are different, 
these structures are separately treated to hi^xLighten the 
important features involved in their design. 

6.1 Pin- jointed structures 

Reference t Chapter V : Illustirative Prohlon No. 1 

This is a plane truss wilh one degree of indeteir- 
minacy and is subjected to two load conditions. It took six 
cycles to achieve the stationary values of the joint displace- 
ments by the point iteration technique with an initial 
assumed zero displacements for all joints. All the member 
areas were assumed to be equal at the outset. The first 
feasible solution was obtained by prorating all the areas of 
the members in accordance with the maximtam member stress. 

This volume was 567*46 cubic inches. The manbers were then 
proportioned to achieve full stress for all members in either 
of the load conditions. For these member areas and joint dis- 
placements, only three eycles of iterations were required^ for 
the point iteration procedure, to converge. The member stresses 
deteamined and member &u?eas prorated to achieve another feas- 
ible solution. Graph 5*1 shows the volume for feasible struc- 
tures versus the cycles of iteration. The minimum volume 
was obtained at the sixth iteration. For this volume only 


r24 

three of the five members were fully-stressed under either 
of the load conditions. The full-stress condition was ob- 
tained at the 15th cycle with a volume of 557*5 cubic inches. 
The difference between the first cycle volume and the minimum 
volume was found to be 5?^ of ihe economic volume. 

Reference: Chapter V : Illustrative Problem No. 2. 

This is a space truss with three degrees of 
indeterminancy and is designed for two load conditions. The 
first feasible volume of the structure was found to be 5^80 
cubic inches. Graph 5*2 shows the volume for feasible struc- 
tures against cycles of design. For this structure also the 
full stress volume was obtained as 2999*55 cubic inches where- 
as the minimum volume was obtained at the second iteration as 
2966,41 cubic inches. At this minimum wei^t only one of the 
members is fully-stressed and all others are stressed lower 
than the permissible stress. The difference between 1st 
cycle volume and the economic volume is 24?^ of the minimum 
volume. 

Reference: Chapter V ; Illustrative Problem No. 3. 

This is a space truss and represents a hemi 
spherical dome in gecmetiy. Since it is a ring structure 
the structure stiffness matrix does not have a band. It 
has a total of 91 members and the degree of statical mde- 
teiminacy is 47* For this structure the feasible areas 
decreased from a value of 21000 to I276O cubic inches 


representing a volume change of 655^. For this structure at 
minim u m weight all the members are fully stressed at either 
of the load condition. 

From the above design examples, the following 
may be concluded. 

(1) At minimum weight all the members of the struc- 
ture need not be fully-stressed. A feasible solution is 
achieved by iterating and prorating the areas. L.S. Schmit (2) 
arrived at the same conclusion by solving the problems using 
mathematical programming. 

(2) Since the stiffness matrix of a truss have always 
strong diagonal elements the rate of convergence to achieve 
the correct solutions is quite fast. For subsequent cycles of 
design where the structure stiffness matrix does not change 
much it requires only two to three cycles for most of the cases 
for a stationary value of joint displacements. This advantage 
cannot be taken of in the direct procedure where the stiffness- 
matrix is inverted to get the joint displacements. So although 
the point iteration may take a little longer time to achieve 
the joint displacements for the first cycle, for subsequent 
cycles, it may be faster than direct inversion. For example, 
for the space truss of figure 5*5 "the total structure stiffness 
matrix will have a dimension of 39 x 39* For joint displace- 
ment this 39 X 39 matrix is to be inverted. point iteration 
only 13; 3x3 matrices are inverted. If for subsequent cycle 


an average of two iterations are necessary the total number 
of 3 X 5 matrices inverted are 26 where as in the direct in- 
version they are two 39 x 39 matrices. 

( 3 ) By this procedure only storage required for the 
stiffness matrix of the structure are two 3 x 3 02 ? 2 x 2 
^oint stiffness matrix for space and plane frame. This saves 
ccanputer memory to a great extent and even a small computer 
like IBM 1620 can handle structures with hundreds cf degrees 
of freedom. This is especially true for problans like 
Illustrative problem 3 of Chapter 5 where there is no band 
formation in the stiffness matrix. 

(4) It is seen that ttie percentage difference between 
the initial design volume and final design volume differs 
according to the degrees of indeterminacy. This is expected 
as member areas are influenced by redundants. 

(5) Proration of areas gives several feasible solu- 
tions and generally takes less number of cycles to achieve 
the minimum weight conditions than necessary to achieve a 
fhll— stress design. Por problem number 1» number of design 
cycles were only 6 and for problem 2 it is only 2 to achieve 
the minimum weight where as the full stress design takes 
fifteen cycles for problem no. 1 and over 3® cycles for 
problem 2. But for structures where fdlly-stress design is 
also the minimum weight design, such conclusions cannot be 
applied. 



(6) For externally determinate structures by going 
through the design of members of the indetenainate complexes 
computer time may be saved. 

Althou^ the illustrative examples are all inde- 
terminate structures the' computer programme developed can 
also handle determinate structures. 

The general com'puter programme developed solves 
both space and plane trusses, "the second as a special case of 
the first. The dimension of the structure is represented by 
KD which takes the value 2 for a plane truss and 5 a space 
truss. The programme is in terms of HD and either a plane or 
a space truss is analysed according to if HD is read as 2 or 
5 . The number of stataaents of the programme for joint dis- 
placements are only 39 plus ten or five according to whether 
it is a space or a plane truss. 
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6.2 Rigid jointed Plane Prames 

Reference » Chapter V : Illustrative Prohlem Mo. 4 

This is a square portal, subjected to a hori- 
zontal load of 10 tons at the floor level. Without the 
introduction of initial jeint deformation and for arhitary 
member proportions the convergence for sway displacement 
was extremely slow although the rotations and axial defor- 
mations converged within first four cycles of iteration. The 
sway defomations of the joints are achieved through the 
axial deformation of the beams. So with an arbitary area of 
the beam the pTOcess of convergence became extremely slow. 

For example the final sway deformation , was built up fran 
zero by P.L/AE the axial deformation of the beam where P. is 
the imbalanced sway force at the ith cycle of iteration and 
these are related by the relation 

.Jc 

A = P.L/AE or (l/AE) Z1 \ •* ^*2 

i=1 

where k is the total cycles needed for stationary values of 
the joint deformation. 

It is well-known that the area of the beams do 
not influence the final sway deformations of joints to any 
appreciable order of magnitude. The sway deformation prima- 
rily depends upon the moment of inertia of column group of 
a storey* So this is a very peculiar process as the 



convergence of sway displacements "but not the final sway dis- 
placement depend upon the area of the teams. So with an 
artitary area of the teams the process of convergence of sway 
deformations may greatly te influenced. 

The initial sway deformations and memter sizes 
are calculated ty approximate methods. When these initial 
guess values for sway deformation and memter proportions are 
assumed then it took very few cycles for the convergence. For 
example the approximate sway deformations was 0.60 inches 
where as the actual sway deformation was O .64 inches. The 
final deformation was found on the hi^er side tecause the 
teams are assumd relatively rigid than the columns in the 
cantilever analysis. 

The design area of memters are determined at 
their economic depths and for a particular set of loads and 
a span there are several areas tut the minimum one was achieved 
ty iteration on depidi. For columns this iteration converged 
within four to five cycles where as for team it takes 5 to 15 
cycles to converge depending primarily upon the tending moments 

moments. 

For this case the colxnnn proportions converged 
with in the second design cycle within an accuracy of 
whereas the teams took two more cycles to converge. 

The sway deformation at toth tte ends of the 
team (namely at joint 1 and 2) remain same except for a 
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small percentage (5?^) due to the axial defoimation of the 
beams* 

Eeference s Chapter V s Illustrative Problem Ho* 5 

This is a continuous gable with four bays and 
is subjected to both wind loads and graviiy loads. 

The joint defomations for this p3X)blem converged 
faster than the portal# This is because it has several in- 
dined members and the off diagonal term 6 EI/l^ in the shear 
equation which caused poor convergence was distributed throu^- 
out the shear equation which improved the convergence of joint 
displacements. 

Like the previous example to achieve the eco- 
nomic depths the beams take more number of cycles than the 
columns. 

In the design cycle all but the first a33d last 
column size converged wi-Hiin three cycles whereas the beams 
took 5 cycles. to converge. This is because both these columns 
carry hi^ moments as the horizontal loads are concentrated 
at these member ends. Both fifth and sixth cycle member 
sizes are given in the illustrative chapiter 5* 

Reference : Chapter V s Illustrative Problem No. 6 

This is a five bay and five storey building 
frame. The process of achieving the final design was com- 
pletely autcanated. Only the number of bays, number of 



storeys and the span of each bay, height of each storey, the 
joint loads and joint coordinates are the input data. Rest 
of the calculations are done within the computer. 

For this problem the sway deformation inti?o- 
duced by approximate method was helpful but the final and 
initial values different by 10C5^ for the top floor as assumed 
sway for top floor was 0.001 5 whereas actual sway was 0.0029 
inches. For bottom floor the respective values are 0.00039 
and 0.00099* These are small values and are due to sway 
and axial deformation of beams. Hence although the results 
for this problem does not look encouraging, since these are 
small values and their contribution to the force equilibinum 
equations are not significant, the convergence was faster. 

For example all the deformations converged within the first 
20 cycles. 

In the design qycle the column convergence was 
foCind to be faster than the convergence of beams. The ini- 
tial design and final design agree for columns within 4C^ 
whereas beam sizes change from 50 to 100^^ between prelimi- 
nary design and first cycle design. 

Prom the above design examples the following may 

be concluded. 

(l ) As alread7 noted for trusses the storage is only 
two 3x3 stiffhess matrices for any "type of plane frames. 

By this method even using a auall computers like 1620, 



structures with hundreds of degrees of fireedon, can he solved. 

(2) Like trasses structures without a proper band 
formation can be solved in this method whereas the solution 
of such types are hi^ly involved by regular procedure. 

(5) The operations . and storage in a computer can 
still be eoononiaed if sway for a storey is assumed constant 
for all the joints of the storey and axial deformation of 
columns are neglected. Columns carry predominantly the 
axial load and if the axial stresses are same in magnitude 
and sense for all the columns of a storey the strain and 
also the axial deflect^^Lons of all the columns of the stor^ 
with uniform column heights will be same. This means the 
beams will sink uniformly creating no beam forces due to 
axial deformation. But in case of wind loads where columns 
are in tension such conclusions oaamot be drawn and axial 
deformations are to be considered. In the design, the ele- 
ments have the absolute minimum weight under the loads. This 
saves the weight to a hi^ degree. The depth of an element 
solely depends upon the forces it has to carry. The depth 
of 1/10 of the span generally assumed, although may be cor- 
rect for (juite a few cases but may be uneconomical for 
several other cases. For example for beam 1 (Fig. ) for 
a depth of I/10 the span the area is 55 square inches where- 
as the minimvim area is 33*62 square inches at a depth of 
26.74 inches. For this case the design office assumption 
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is corrects But for beam 2 the minimum area is I6.4 sq,. 
inches at a depth of 15.4"* If a depth of l/lOth the span 
is assumed, the area is as hi^ as 19«58 sq. inches. So 
in this case the increase in area is nearly 269S. So if tl| 
areas are calculated at the economic depths 20?^ of the toti 
frame wei^t may be saved. 

Once the section areas and mcOTent of inertia 
are thus obtained the sections can either be fabricated or 
be selected from handbooks. Economic factors like fabrioa 
tion should govern the selection. 

This way of design of a structures has two 
distinct advantages. I 

(1) It saves material as the elements are designed 
for minimum area. 

(2) The structure designed is safe and at no place ’ 
will be overstressed. In the usual design where a structure 
is analysed and the maabers sxe proportioned for the analysis 
the final design stresses may not be same as the assumed 
stresses an^ hence in some cases overstress of the manber 
may result. 

6.5 Future Work on this Field 

The structures designed are elastic. But as 
load increases some of the mehers oome to the plastic zone. 
So Instead of going on an elastic path, the stress-strain 
curve may be asstjmed non-linear and at the membOT load 
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proper elastic modulus ’E’ may be taken according to the 
actual case. This will change only the STJBSOTJTIKE SDIFIJ 
and STIFJ- which calculate the elem^t stiffness matrix in 
the local coordinate system. 

Once the member sizes are detemined, these 
members may either be fabricated or for the areas and 
moment of inertia standard sizes may be taken from proper 
handbooks. The section which gives maximum economy may be 
the governing criterion. For example if the fabrication 
cost is higher these standard sizes should be used otheair- 
wise they may be fabricated. 



CHAPTER VII 


EXPLANATION OP COMPUTER PROGRAMMES AND FLOW CHARTS 
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M&IN FBOGEAMME 


mPDT 

(i) Geometry of the Frame 

The total number of joints of the frame; the 
number free joints; Number of Bays and Floors; Hei^t of 
each Storey and Span of each Bay; Coordinates for all joints. 


(ii) External Loads 

Forces acting at each joint, load and type of 
load along the span of each member. 


EXECUTION 

CALL APPROX.* 


Analyse the structure by stiffness method. For 
each joint the joint stiffhess matrix is generated. The si mu l' 
taneous eq,uations are solved by point iteration procedure and 
joint displacements are determined. The above procedure is 
repeated for all joints. Structure equilibrium is checked 
at each cycle snd iteration is said to be stationary when the 
structure equilibrium is satisfied. Member end forces are 
calculated and design forces determined by CALLING FUNCTION 
COUPLE, SFORCB and AXIAL. DESIGN SUBROUTINE is CSJitS® and 
all membeirs are designed. The structure is reanalysed for 
these member sizes. This iteration is repeated until the 


member sizes ranain constants at two consecutive cycles. 

♦"aPROX is a^UBRODTINEV tte function of it is summarized 
under the SUBROUTIHE Catalogue# 



CfUTFOT 


Depth, Areas, Maaent of inertias for each member 
and the value of interaction fomula. 


CATALOGFDE OP SUBEOUTINES 

SUB 1 1 SUBROUTINES A P R 0 X 
INPUT 

The number of Storeys and Storey heights. 

Number of Bays and Span of each Bay. 

Number of joints and coordinate for each joint. 

Sway forces at each storey level. 

Load and Load type on the Span of each Beam and 

Column. 

EXECUTION 

Calculates the number of manbers meeting at a joint 
and the far end joint numbers for each joint. 
Calculates Hie Bending Moment, Shear Force, Axial 
Force for each member due to wind forces assuming 
the building as a Cantilever and the same due to 
gravity loads assuming the beams to be pinned to 
the Columns at both the ends for Column loads and 
completely fixed at both the ends for beam loads. 
Adds up both the effects and determines the maaber 



end moments y sheax forces and axial forces. The 
design forces for each msnber is obtained by 
using Bending Moment (FDMCTION GOUHiE) , Shear 
Force (FOTCTIQN SFORCE) and Axial Force (POMC- 
TION AXIAL) diagrams. 

Each member is proportioned for minimum area 
by CALLING the DESIGN SDBBOTJTINE. 

For the above member forces and sizes, the sway 
deformations are calculated (PONCTION DELTA) at 
each floor level assuming floors to be rigid 
and the appropriate joints are assigned these 
values. 

ODTPDT 

The number of members meeting at a joint and the 
far end joint numbera for each joint? the areas, 
moment of inertia for each member and the sway 
deformations for each joint. 

SOB 2* SOBROOTINEi DESIGN 
INPDT 

Design Mcanent, Shear and Axial Force; Span; 
and permissible Stresses in bendiaig, axial 
and shear. 

EXECUTION 

Assumes a depth 1/20 of the Span to satisfy 
serviceability. For th^ depth designs an I- 
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criteria. Changes the depth by ^ and repeats 
the above cycle until the economic depth is 
reached, i.e. the depth at which the areas starts 
increasing. 

OUTPUT 

Member area, Moment of inertia, Depth and the 
Interaction Value for the section. 

SUB 5* SUBRCfUTIMEs M A T M P 
INPUT 

Matrix A and Matrix B 


EXECUTION 

OUTPUT 

SUB 4* 
INPUT 


Multiplies A and B and stores in matrix G 


Matrix C 

SUBEOUTINEs MARTIN 


DENCM, a, 3 X 3, Matrix 


EXECUTION 


OUTPUT 


Inverts DENCM and Stores it in BINV Matrix 


Matrix DINV 

SUB 5! SUBROUTINE: TRANS 

INPUT 


The Direction Cosine of aaiy member. 



EXEODTION 



Generates the IMSFOIMATiaEI Matrix ST for the 
member. 


OTJTPOT 

Transfomation Matrix ST. 

SUB 6t SDBROUTIKE: X L 0 A I) 

INPDT 

Joint Loads, Span Load, Joint Coordinates, 
Total number of members meeting at any joint 
and the far end joint numbers for each joint. 

EXECUTION 

Determines the equivalent joint loads due to 
inspan load and then adds this to the loads 
directly applied to the joints to give the 
Joint Load ?eotor. 


OUTPUT 

Equivalent JOINT LOAD VECTOR 
SUB 7* SUBRDUTIHEt S T I P I J 
INPUT 

Member Areas, Moment of inertias. Span and 
Young’s Modulus of the manber material, near 
and and far end joint nximbers of the member. 

EXECUTION 

Generates SK(KjP matrix in the local coordinate. 

OUTPUT 


SE Matrix 
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SUB 81 SUBEDUTINEs S T I P J I 
IHHTT 

Same as SOBEOUTIBE STIPIJ 

EXECUTION 


Generates SSK(K^j^) matrix in the local coordinates 


OUTPUT 


SSK Matrix 


SUB 9 s FUNCTION; DELTA 
INPUT 

Bending Moment » Span and moment of inertia of 
members > gene3?ally the Columns and Young’s 
modulus of the material. 

EXECUTION 

Determines the individual SWAY Deformation 
as 6 1 , 1 / 1 ? 

OUTPUT 

Sway Deformation 

SUB 10 s FUNCTION; COUPLE 
INPUT 

Member End Moments, Shear Forces, and Axial 
Forces. Type of Load on Span, and distance, x, 
to a section from the near ©ad A for which all 
forces axe known. 


EXECUTION 

Determines the bending maaent at ary section 
within the Span. 
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OXITPDT 

Bending moment at any section along the span. 

SOB 11: FONCTION: AXIAL 

INPUT 

Same as .FUNCTION COUPLE 

EXECUTION 

Deteimines the axial force at any section 
within the span. 

OUTPUT 

Axial Force at any Section. 

SUB 12: FUNCTION: S F 0 R C E 

INPUT 

Same as FUNCTION COUPLE 

EXECUTION 

Deteimines the Shear Force at any section 
along the span. 

OUTPUT 


Shear Foixse at any section. 
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CCMPDTER PEOGEAMME POE TRUSSES 

INPUT 

Total number of joints of the frame. Dimension 
of the frame; Number of free joints; Number of 
Boundary joints and degrees of freedom at each 
Boundary joint; Joint Coordinates; Joint Loads; 
Par end joint number for each joint. 


EXECUTION 

For all monbers meeting at a joint generates 
the span) the direction cosine and transformation 

m 

aatrix ’T’and the product (T T)# Member stiff- 

T 

ness is determined by multiplying the T T 
matrix by AE/L. 

The process is repeated and joint equilibrium 
equations are generated using the above member 
stiffness matrix according to the relation 2.11 
and joint displacements calculated according to 
equation 2.12. 

The process is repeated for all joints several 
times till the joint displacements have a station- 
ary value. 

Member end forces determined according to equa- 
tion 2.14. The above procedure is repeated for 


all load conditions. 
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For multiple load condition memlDer areas are 
prorated according to the maximum mmber stress 
if this is more than the permissible value 
otherwise members are proportioned to achieve 
full stress. 

For single load condition repeat the above pro- 
cedure accept that only (M-N) member are fully 
stressed. 

OUTPUT 

Member area* and member stresses » volume of 
structure for each design cycle. 
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FLOW CHART FOR THE LESIGH OF A. RIGID JOINTED 
FLARE FRAME 


/ 


START 


A 


M, R, IR, RD, IFLOOR, IBAY, FORCE, P, W, X, 


COLDl'lH, DL, ASPAR, SPREAD 




’ScritializatiohI 


, DO = 0 


i GENER/iTE JTOT ARD JOIRT MATRICES ! 


FOR LATERAL LOAD ARALYSE THE STRUCTURE BY 
CARTILEYER METHOD ARD FOR GRAYITY LOAD 
ARALYSE IT ASSDMIRG BEAMS SIMPLY SUPPORTED 
FOR COLUMR REACTIOR ARD RIGIDLY FIXED FOR 
BEAM MOMENTS? DETERMINE DESI® LOAD 


'CiLrDESiOR’SUBEOUTIHE ' 
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I 

} 

{ 

POE THE MEMBEE SIZES; DETEEtUHE : 

I 

4 PPE 0 XIMATE SWAY DEPOBMaTION ; 

I 



AHILYSE THE STEUCTHEE BY POIM 
HEEATIOE OP THE STIPFlffiSS 
EQUATIONS 



v’ * 


INITIALIZATION 

INITIALIZATION 

DENOM = 0 



UPPEE = 0 



GENEEATE SPAN, TT, ST, SST Jlatrioes 

DENOM = XSST ST; UPPEE =^SST ST D(J) 

I WEST, WQfl . IN . sml 

I" D(l)'*'=’'Dnw"» (P(l)-UPPPl) I 


0 
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DETEacHJE TOTAL UNBALAHCED 
JOIMT PQECE AKD CHECK 
TOLEEAITCE 



0— '■ 


:uo = u 

I 

— ) ^ j - 

Ciir'cOuiLE*, SPORCE, mAL; 

— t 





-YES 


1 


YES 

I 

f 


;' PRTfIT SECTION 

: iifitfiEfes ' 



STOP \ 

i 

y 
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COMPUTER PROGEMME FOR PIN CONNEX3TED FRAMES 


/ N 

START ' 


M, N, IN, ND, JOIM’ ' (“I ,jyr JTOrri) , 
X, P, IPROB 

r 


• INITIiLIZATION . 

I 

U « 

0 ! 

I uo « 

0 I 

f” ‘ ‘ 

IL = 

1, IPROB 



I =s 

“fn 

1 - 1- 1 -1 r 

, _ , , 




J =1, jtot(i)! 

“ki' 'i j6iHT(i,J)' 

GENERATE T AEB T^. I 

PRODUCT T ■* T^; DENOM = DENOM + J 
DENOM = DENOM + T^ * T^ * AE/L ! 

» UPPER + T * T^*AE/L j 

( 
i 

J 


UPPER 
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d(i) = denom”"’ (p(i) - iipper) 



I MEMBER EORCES ACCOEDUG TO EOa’IULA 



D^EMHE JlBffiER ST^SS 



NO 



YES 


O'/ERSTEESSED 




PEOP.ATE AREiS 


, PROPOECION 
; AREAS TO ' 
ACHIEVE i 
. POLL STRESS I 


PRIHT MHv'TBEE FORCE, AREA, STRUCTURE i 
VOLuilE ■ 
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NOTATIONS USED IN CCMFDTER PBOGMIIMES 


A, AAREA, ABEA 

AIXX 

ALFA 

ASPM 

AX 

AXMAS 

BBB 

BEAM 

BEAM 

BEND 

BETA 

BEL 

MA 

BMB 

BRAT 

CAIXX 

CiJREA 

CG 

CHS 

CO 


Area of the members. 

Moment of ineirtia of members. 

Tolerance matrix for load 
Span of tte Bays 
Axial "FoToe 

Design value of axial force 
Intemediate matrix used to convert local 
in span load to global joint load 
Moment of beam due to wind load. 

Shear force of Beam due to wind load 
Fraction of bending stress used in the 
interaction formula 
The interaction ratio. 
ld/b«t 

Bending moment at end A 
Bending moment at end B 
Ratio (Ld/bt)/600 
Column moment of inertia. 

Column area 

Centre of gravity of the column group 

of a storey. 

Column Tinit shear force 
Bending moment at any section of 


a member^ 



COLD 

COLP 

COLM 

COLS 

COLDMH 

COLTO 

COMAX 

COOPLE 

DEI/PA 

DEN 

DENO 

DENCM’ 

DEPTH 

DIP, DIPPE 

DIN7 

DL 

DP 


E, Ei 
ERA, EEB 
F2 


CoxuQui axial force due to ■<=^avity 3oad. 
Column axial force due to gravity load. 

“ Column moment. 

Col\mm shear force. 

“ Height of columns. 

Column normal force. 

~ Design moment for any member. 

“ Function which draws beaiding moment 
diagram. 

Sway displacement. 

Intermediate matrices in the generation 
of total denominator 

Depth of section. 

- Difference of joint displacement between 
two consecutive cycles. 

“ Intwrse of joint stiffness matrix K. . 

Dead load on span. 

- Maximum depth to prevent holding budding 
of the flanges of a beam. 

- Youngs modulus of the material. 

- Axial force at aid A and B of a member. 

- Intermediate matrices for evaluating 


member end forces 
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PA2 


Flange area for both flexure and axial 

force. 

PAL 

P. 

The axial stress to prevent compression 

buckling. 

PM, PEB 

- 

Fixed end moments for end A and B. 

PEM 

- 

Fixed end moment. 

PLAUG 

- 

Design area of flange. 

PMAX 

- 

Design thickness of flange. 

PORCE 

- 

Sway force at storey height. 

IWIDTH 

- 

Design width of flange. 

IBAY 

- 

Number of Bays. 

IBL 

- 

Beam number at left hand side. 

IBR 

- 

Beam number at rigjit hand side. 

IBEAM 

- 

Total number of beams. 

ICOLUM 

- 

Total number of columns. 

ICYCLE 

- 

Iteration cycle for design. 

IDOWTIL 

- 

Left hand side joint number for storey 

below. 

IBOWNR 

- 

Right hand side joint number for storey 

below. 

IPIBLB 

- 

Total number of columns at ai^ storey. 

IPBOOR 

- 

Total number of floors. 

IN 

- 

Maximum number ff members meeting at 

any joint. 

laERXL 

- 

Left hand side number of the joint 


of a storey 



ISTETR 


Ri^t hand side number of the joint of 

a stress. 

ITERA 

- 

Iteration cycle for analysis. 

lOPL 


Left hand side joint number for storey 

above. 

lUPR 


Right hand side joint number for storey 

below. 

JOINT 

- 

Par end joints for any joint. 

JTOT 


Total number of members meeting at any 

joint. 

KI 

- 

One far end joint. 

M 

mm 

Total number of members including boimdary 

joints. 

M> 

- 

Top number of a column. 

m 

- 

Bottom number of a column. 

N 

- 

Total free joints of the structure. 

ND 

- 

Dimension of structure. 

P 

- 

Loads at joints. 

PP, PPR 

— 

Fixed edge forces in the member In- 

system coordinate. 

RMIN 

- 

Minimum radius of gyration. 

RXX 

- 

Radius of gyration along X axis. 

RTY 

- 

Radius of gyration along Y-axis. 

S 


Difference of joint coordinates of 


a member 



SALOT 

SF 

SFOEGE 

SH 

SHEAE 

SHMAX 

SIGMA 

SK 

SPAN 

SPREAD 

SSK 

ST 

STT 

SWAY 

T 

TP 

TEMAX 

TOLBNC 

TOTLOD 

TT 

TW 


•" Allowable shear stress in the web of a 
beam. 

“ Unit shear force. 

- Function which gives the shear force 
diagram. 

“ Shear force at any point of a nmember. 

- Shear force at any storey level. 

Design shear force. 

Pemissible stress in bending or axial 
compression. 

~ K. .-matrix in local system. 

i J 

Span of member. 

- Span of beams. 

matrix in local coordinate system. 
Systoa transfer matrix. 

- Transpose of transfer matrix. 

Sway displacements. 

Direction cosines of a manber. 
Thickness of flange. 

Design thickness of flange. 

- Allowed tolerance in the equilibrium 
equations. 

Total load acting in the structure. 

- Direction cosine of ai^r member. 

- Thickness of web. 



TTOlftX 


^56 


V 

UK, UKK 

UKBL 

UO 

UPP 

UPPE 

UPPER 

UPPT 

UU 

VALOW 

IW 

X 


- Design thickness of web. 

- Joint displacements. 

- Matrix used to generate the numerator 
of the joint stiffness matrix. 

- Total unbalanced load at any cycle. 
Initial displacaaents. 

I 

^ - Same as UK? TIKK 

11 « 

j 

- Allowable shear force in the web of a 
beam. 

- Span load. 

- Joint coordinates. 

Moment of inertia of a member. 


ZI, ZXX 
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